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ABSTRACT 


This thesis utilizes Military Specification, Mil-A-8591H, 
to provide and document a new computer program for calculating 
store reaction forces, due to inertial loads and airloads, at 
the aircraft-store interface. Due to the fact that the store 
is a statically indeterminate structure with unknown 
flexibility characteristics, the analysis to calculate the 
reaction forces is a relatively complex problem. The program 
was written to incorporate new requirements and techniques for 
calculating reaction loads at the aircraft-store interface in 
accordance with Mil-A-8591 revision H. This computer program 
was developed to be used on the VAX/VMS Computer System, 
Advanced Computational Laboratory, of the Aeronautics and 
Astronautics Department, Naval Postgraduate School. A detailed 
User’s Manual which is included in Appendix A of this thesis 
was written to accompany the computer program. The computer 
program will enable graduate students in the Aircraft, 
Helicopter and Missile Design courses to utilize current 
military specifications for preliminary design and weapon 


system integration. 
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THESIS DISCLAIMER 


The computer program waS written to provide a 
Simplified analysis for preliminary design. The program was 
Written according to the methods and equations provided by 
military specification but has not been tested for all cases. 
Any application @is tars program without additional 


verification Sseatrene risk of sene User. 
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I. INTRODUCTION 


The design of a new military aircraft is predicated on 
wnat stores the aircraft is required to carry. A major design 
consideration for aircraft and aircraft systems is the 
interface between the aircraft and the store. 

Due to the fact that the store is a statically 
indeterminate structure with unknown fitewubi lity 
Characteristics, the analysis to calculate the reaction on the 
store suspension equipment is a relatively complez problem. By 
making several assumptions, the problem can be simplified so 
a general analysis can be performed for stores and aircraft 
that are still in the preliminary design stage. 

Military specification MIL-A-8591H is the general guide 
fem the integration of a store to an dircraft. It covers 
requirements that need to be considered in a store-aircraft 
Saeegration, including the main focus of the specification of 
how inertial loads and airloads affect the aircraft store 
maeterface. Mil-A-8591H, Appendix D (Ref. 1], outlines an 
updated approach to analyze the results of specified inertial 
loads on a store-aircraft suspension equipment interface. The 
objective of this thesis is to provide a computer program and 
ieer’s manual for calculating the reaction at a store-aircraft 


interface due to specified inertial loads and airloads in 


accordance with the method described in Mil-A-8591H, Appemcam 
D [Ref. 1]. 


II. MILITARY SPECIFICATION MIL-A-8591H 


Mil-A-8591 was first published in 1955, since then it has 
gone through several revisions until the most current version, 
Mil-A-8591H was published on 23 March 1990. One of the primary 
goals of the specification is to provide aircraft store 
standardization not only for the United States Department of 
Defense but also for NATO (North Atlantic Treaty Organization) 
countries. The specification contains detailed drawings for 
store hook and sway brace design requirements. The 
specification also provides general guidance for design, 
analysis, test and documentation of airborne stores, 
suspension equipment and the aircraft-store interface during 
captive operations. In addition to the general body of the 
specification there are four appendices that give detailed 
guidance to analyze the effects of inertial loads and airloads 
on the aircraft-store interface. 

Mil-A-8591H, Appendix D [Ref. 1], provides an indepth 
method to analyze the reaction forces at the aircraft-Sstore 
interface. Through the years there have been several analyses 
for the reaction forces at the aircraft store-interface. One 
of the most recent is contained in [Ref. 2]. The method of 
analysis in Mil-A-8591H was not addressed in previous 
revisions of the specification. The various analyses rely on 


basically the same information (1.e., store physical 


characteristics, store inertial loads, and store airloads) to 
provide reaction forces at the aircraft-store interface. The 
analyses are conservative and cover a broad category of stores 
and a wide range of aircraft. Once design considerations 
become more specific a more detailed analysis can be 
accomplished. 

Mil-A-8591H, Appendix A [Ref. 1], provides store inertial 
loads and store airloads for the most general analyses. The 
store inertial loads for fixed wing high performance carrier 
based aircraft are provided. Also given are equations that 
provide maximum store angle-of-attack and maximum store angle- 
of-sideslip so that corresponding airloads can be calculated. 
As would be expected, the approach in Mil-A-8591H, Appendix A 
[Ref. 1] is the most conservative and covers a broad band of 
aircraft and stores. Note, that in Mil-A-8591H, Appendix A 
[Ref. 1], store envelopes are provided, yet no reference is 
made to aircraft envelopes. 

Mil-A-8591H, Appendix B [Ref. 1], 1S more specific in that 
the inertial envelopes of specific types of aircraft are 
provided. In this case the store inertial loads must be 
calculated using the accompanying aircraft equations of 
motion. The analysis in Mil-A-8591H, Appendix B [Ref. 1], 
requires more specific information that is not always 
available during preliminary design. Since the analysis of 
Mil-A-8591H, Appendix B [Ref. 1], is for a specific aircraft- 


store combination, the results are less conservative than the 


results of the analysis of Mil-A-8591H, Appendix A [Ref. 1]. 
Meme, in Mil=A-@591H, Apoendix Bai Rets 1), in contrast to Mil- 
A-8591H, Appendix A [Ref. 1], aircraft envelopes are provided 
and the store envelopes must be calculated. 

Mil-A-8591H, Appendix C [Ref. 1], is specifically for 
helicopters and the aircraft inertial envelopes are provided. 
The necessary equations to convert the aircraft inertial loads 
to store inertial loads are in Mil-A-8591H, Appendix B 
thet. 1). Specific detailed information about aircraft and 
Store physical characteristics are required. Mil-A-8591H, 
Appendix C [Ref. 1], provides no specific guidance to 
calculate airloads for a helicopter-store combination. 

All analyses and envelopes provided in the specification 
are general in nature and all stores shall conform to the 
requirements of the specification, but the acquiring activity 
could prescribe more stringent or detailed requirements in the 


aircraft or store detailed specification. 


III. AIRCRAFT STORES INTERFACE MANUALS (ASIM) 


The ASIM manuals are a reference to provide engineers a 
source of information regarding the phySical characteristics 
required in aircraft-store integration. The information 
provided by these manuals is essential for new aircraft 
design, new store design, and existing aircraft-store 
integration. Since future aircraft and aircraft system design 
builds on what already exists, the ASIM manuals are a good 
Staverng: perme for any  aircraft-store analysis. The 
information in the ASIM manuals will be utilized in the 
Mil-A-8591H loads analysis. The ASIM manuals are divided into 
three volumes [Refs. 3, 4, and 5]. See Figure 1 for the 
Organization of ASIM. 

Volume 1 is the Aircraft Manual and it provides detailed 
drawings of current United States and NATO aircraft store 
stations. The manual is arranged into four sections and an 
inGsoductery SeGE lon: 

The Introductory Section of Volume 1 ties together the 
other four sections by describing the contents and suggesting 
typical uses of the manual with an overview of store 
integration. In addition a chart is provide to indicate tae 
SuspenSion equipment currently being utilized on United States 


alrcrart. 


Aircraft Stores Interface Manuals 
(ASIM) 


Volume 2 Volume 3 
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Figure 1. ASIM Organization 


Section 1 of Volume 1 of the ASIM manual provides bons 
rack drawings with associated dimensions for integration to 
the aircraft. Similar information is provided in Volume 3, 
only in much more detail and with respect to the entire bomb 
rack or suspension equipment. Section 2 of Volume 1 of the 
ASIM manual provides aircraft drawings of United States and 
NATO aircraft. The aircraft are presented giving a front view 
and side view together with dimensions required for store 
integration. Section 3 of volume 1 of the ASIM manual provides 
composite drawings of several aircraft so that clearance 
envelopes can be compared. The composite nature of the 
drawings allow a designer to determine which aircraft can 
accommodate a store and what clearance iS available for a 
Store designed for use on a variety of aircraft. Section 4 of 
Volume 1 of the ASIM manuals provides detailed drawings of the 
aircraft pylons to which the stores must attach. 

Volume 2 is the Store Characteristics Manual and it 
provides detailed drawings of current stores utilized by the 
United States and NATO. The manual is broken down into 14 
sections. The sections cover, Conventional Bombs, Nuclear 
Bombs, Practice Bombs, Flare Dispensers, Rockets, Air-To- 
Ground Missiles, Air-To-Air Missiles, Mines, Torpedoes, Guns, 
Fuel Tanks, Pods, and Targets. The primary dimensions, weights 
and mass moments of inertia are given for each store. In 


addition, the phone numbers and addresses of the design and 


field activities for each store are listed, if further 
information about a store is needed. 

Volume 3 is the Suspension Equipment Manual, and it 
provides detailed drawings of the suspension equipment used by 
United States and NATO countries. The manual is broken down 
into sections describing various types of bomb racks and 
missile launchers. A detailed description of the suspension 
equipment capabilities and performance characteristics is also 
included. Also given are phone numbers and addresses of the 
design and field activities for each store. These are listed 


if further information about a store is Sought. 


IV. PROBLEM FORMULATION 


A computer program has been developed to implement the 
store interface analysis contained in Mil-A-8591H, Appendix D. 
The analysis of Mil-A-8591H, Appendix D [Ref. 1], was not 
published in previous revisions of Mil-A-8591. Existing 
programs [Ref. 2] use a different approach in the analysis. 
Reference 2 analyses tandem hook/lug tandem sway brace 
configurations and single  hook/lug (practice bombs) 
configurations. The basic premise of the following analysis is 
a tandem hook/lug tandem sway brace configuration, Figure 2. 
The objective of the analysis is, given a specified flight 


envelope, to determine the worst case load at the aircraft- 


store interface (i.e., hook/lug and sway braces). 

The basis of the analysis is a set of equations to 
determine the worst case load for a component given the 
inertial loads as specified in Mil-A-8591H. If airloads are 
applicable they must be converted to forces and moments about 
the store CG and added to the corresponding inertial loads. 
The computer program allows direct Input of store 
accelerations due to aerodynamic effects and has the 
capability to perform a basic airload analysis from References 
Lj ib sanahas 

The overall problem is complex with many variables not 


available during preliminary design. In order to reduce the 
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Figure 2. Relationship of Typical Bomb-Rack Top-View 
Drawing from Section 1 to the Actual Rack and the 
Carrying Aircraft 


‘1: 


complexity of the problem (i.e., to a statically determinate 


structure) several assumptions are postulated. These are: 


Treat the store as a rigid body. (Assuming the store a 
rigid body, the flexibility of the store, resultant 
load paths and modal bending are removed from the 
analysis.) 


A sway brace can only react compressive loads, while a 
hook/lug can react both tensile vertical loads and 
Shear loads. 


The designer is only concerned with worst case loads on 
the store attachment points. (By analyzing each sway 
brace and hook/lug, utilizing the sign convention of 
Figure 3, the worst case store inertial loads for each 
store attachment component are listed. See Figure 4. 
For example, for a left front sway brace, the worst 
case loading is when the store is accelerating forward 
(-NX), accelerating to the left (-NY), accelerating 
upward (+NZ), roll accelerating clockwise looking 
forward (-PHIDD), pitch accelerating nose up (+THEDD), 
and yaw accelerating nose left (+PSIDD) .) 


Worst case inertial loads are to be applied at the 
center of gravity, CG, of the store. The dimensions of 
the CG to the store attachment points are defined in 
Figure 5 and are described in Appendix A. From this a 
Summation of the forces, 2 Fx=0, 2% Fy=0, 2 Fz=0 and 
Summation of the moments £ Mx=0, 2% My=0, 2% Mz=0 can be 
calculated. The results are the worst case loads at 
each component as a function of the store inertial 
accelerations. 


The equations for the worst case load at the FORWARD SWAY 


BRACE are: 


Vertical forward sway brace reaction; 


SZEWDY=SA*PX + B*PY + C*PZ + (Ose eis 


where: 


PX, 


PY, and PZ are the total forces (inertial and 


aerodynamic) in the X,Y, and Z directions respectively. 


FZ 
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Load-Factor Envelope for 
a p Wing Mounted Stores 
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Figure 3. Coordinate System, Sign Convention, and a Typical 


Load Factor Envelope 
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NOTE: (1) If the CG of the store is located laterally outside 
of the SB’s, a negative vertical load (P,, down) may be 
Critical and should be investigated. 


Figure 4. Direction of External Loads and Moments’ for 
Maximum Reaction Forces at the SB’s/Stores SB Pads 
and Hooks/Lugs 
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Figure 5. Store Physical Characteristics 


Js 


MX, MY, and MZ are the total moments (inertial and 
aerodynamic) about the X, Y, and Z axis respectively. 


Lateral forward sway brace reaction; 
SYFWD = SZFWD * TAN(EPS) 
Total forward sway brace reaction; 


STOTFWD = SZFWD / COS(EPS) 


where: 
Lh = GS m js} 
2(Sl+L2) S* TAN(EPS) 
BR SZ * 15 
S*K 
C = | Se} *(.5 +] 
SS K 
D _ SZ 
S*K 
E = a) 
Su ee 
F aa 
TAN (EPS) 
H = (£5 —So.)) * TAN CEES) 
K = H+L3+S3 


note: SZFWD? + SYFWD? = STOTFWD? 
The worst case load for the left or right forward sway 


brace is the same due to the symmetry of the inertial 
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envelopes and the fact that the stores are carried on the left 
@e right wing. 
The equations for the worst case load at the AFT SWAY 
BRACE are: 
Vertical aft sway brace reaction; 
SAG Geek tS BX ts ere + Ds tEAM yt PAZ 
where: 


Eee FY, ana PZ are the total Eorces (inertial and 
aerodynamic) in the X,Y, and Z directions respectively. 


MX, MY, and MZ are the total moments (inertial and 
aerodynamic) about the X, Y, and Z axis respectively. 


Lateral aft sway brace reaction; 
SYAFT = SZAFT * TAN (EPS) 
Total aft sway brace reaction; 


STOTAFT = SZAFT / COS(EPS) 


where: 
a L5 ‘ L3 
2(S2+L1) S* TAN(EPS) 
B - Sl] * 55 
S*K 
Cc = {Se} *[ 5+] 
ne & 
D _ S1 
S*K 
= “A 540 
Seo 


‘lia 


it 


f= SFAN(EPS) 
H = (235-85) * TAN(EPS) 
K = H+L3+S3 


note: SZAFT? + SYAFT? = STOTAFT’ 

The worst case load for the left or right aft sway brace 
is the same due to the symmetry of the inertial envelopes and 
the fact that the stores are carried on the left or right 
wing. 

The equations for the worst case load at the FORWARD 
HOOK/LUG are: 

Vertical forward hook/lug reaction; 

LZFWD = AZ*PX + BZ*PY + CZ*PZ + DZ*MX + EZ*MY + FZ 

Lateral forward hook/lug reaction; 

LYFWD = BY*PY + CY*PZ + DY*MX 

Longitudinal forward hook/lug reaction; 

LXFWD = PX 
Total forward hook/lug shear reaction; 
LTOTFWD = ( LYFWD? + LXFWD* )?” 
where: 


PX, PY, and PZ are the total forces (inertial and 
aerodynamic) in the X,Y, and Z directions respectively. 


MX, MY, and MZ are the total moments (inertial and 
aerodynamic) about the X, Y, and Z axis respectively. 


where: 


gg. L3 


AZ = 
(SZ2+I1) S* TAN(EPS) 


abe: 


BZ 
S¥K 
ee = = JI ee) +1) 
L (H+S4-L3) 
DZ = Je 
S¥R 
EZ = : 
ae 
FZ = : 
“S * TAN(EPS) 
AY = EY = FY = 0 
BY = [52 2 tare) 1] 
S K 
a L2 * L3 * TAN(EPS) 
L* (H+ S4-L3) 
- S2 * TAN (EPS) 
[ — 
H = (L5-S5) * TAN(EPS) 
ioe f+ iS S3 


The equations for the worst case load at the AFT HOOK/LUG 
ere: 
Vertical aft hook/lug reaction; 


ee eae KR + BLAEY + CZ*PZ + DZAMX st VEZ *MY + FZ*MZ 


Wo 


Lateral aft hook/lug reaction; 
EYARE = BY*PY + ° CY*22Z wy 
Longitudinal aft hook/lug reaction; 
LXAFT = PX 
Total aft hook/lug shear reaction; 
LTOTAFT = ( LYFWD? + LXFWD? )1/ 
where: 


PX, PY, and PZ are the total forces (inertial and 
aerodynamic) in the X,Y, and Z directions respectively. 


MX, MY, and MZ are the total moments (inertial and 
aerodynamic) about the X, Y, and Z axis respectively. 


where: 
i en Mees 
(SI*io) S@= TAN (EDS) 
BZ Si. xo 
S¥ kK 
CZ = ce L3 a4 
ie (H+S4-L3) 
DZ = oe 
Sa K 
G7 
Sime) 
me = i 
S* TAN(EPS) 
jo Sega = NG 
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st L5*TAN(EPS) _, 
S = 


BY = = 
K 


L1 * L3 * TAN(EPS) 


CY = 
L* (Hees4 io) 
py = 91 * TAN(EPS) 


S*K 


Bor all theweases listed, thercG of the store is located 
between the sway braces laterally and between the hook/lugs 
longitudinally. If a store has a CG that does not satisfy 
either condition, alternate equations must be utilized for the 
coefficients C, CY, and CZ, see [Ref. 1]. The computer program 
incorporates the alternate equations in the calculation of the 
reaction forces automatically, when required. 

mine direction of the vertical acceleration (PZ) for a 
Critical load on the sway braces can be positive or negative 


when the store CG is laterally outside the sway braces. 


Therefore, both positive and negative vertical accelerations 
must be considered when determining the critical sway brace 
load with the store CG laterally outside the sway braces. 


To incorporate store aerodynamic loads into the analysis, 
the respective aerodynamic loads must be reduced to forces and 
moments at the store CG. If data is available from prior 
sources (1.e., from wind tunnel tests), the data can be 
entered directly. Most preliminary designs are not afforded 


this luxury and must estimate aerodynamic loads. Mil-A-8591H 


ZA 


[Ref. 1] provides a method to determine the store’s maximum 
angle-of-attack and sideslip but does not provide any guidance 
on various aerodynamic coefficients. If coefficients are 
available, the analysis can be completed and the results can 
then be combined with the store’s inertial loads. Hoerner 
(Refs. 6, 7] generated several generic aerodynamic analyses 
for different stores mounted on an aircraft. General values 
for the coefficient of lift, coefficient of side-fones 
coefficient of pitching moment and coefficient of yawing 
moment are provided and utilized in a current store loads 
analysis [Ref. 2]. 

The coefficients that are utilized in the Hoerner analysis 
resolve the aerodynamic forces and moments at the store CG. 
The maximum angle-of-attack and the maximum angle-of-sideslip 
are provided by Mil-A-8591H ([Ref. 1]. Utitlazingaeg. 
aerodynamic analysis of [Refs. 2, 6, and 7] the store 
reference area is computed as the circular area consistent 
with the store radius. The reference length is entered 
directly. If the reference length is not available, a general 
length of 22 * store radius is recommended [Ref. 6]. Since 
most stores are symmetrical from a front view, Y-Z plane, the 
coefficient of lift Slope is assumed to be equal to the 
coefficient of .side-force slope, and the pitch moment 
coefficient slope is assumed to be equal to the yaw moment 
coefficient slope. Hoerner suggests values of .32 for the 


coefficient of lift slope and coefficient of side-force slope; 
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a value of .13 for the pitching and yawing moment coefficient; 


and a value of .3 for the coefficient of drag. The equations 


utilized in the computer program to calculate the aerodynamic 


forces and moments at the store CG are: 


Aerodynamic forces in the respective direction; 


PXA 


PYA 


PZA 


CD * QUES R 
CG * ANGS 2 30UE “Wor 


Cl “ASANGASASCUE BAe oR 


Aerodynamic moments about the respective axis; 


QUE 


ANGA 


ANGS 


0 
C2 ANCA OUE-* SRe~ GER 


C7 ANGSee, CUB * SK * ER 


Lift coefficient slope (per degree) 

Pitch moment coefficient (per degree) 
Side force coefficient slope (per degree) 
Yaw moment coefficient slope (per degree) 
Coefficient of drag 

Dynamic pressure (lb/sq ft) 

Store reference area (sq ft) 

Store reference length (inches) 

store angle-of-attack (degrees) 


Store angle-of-sideslip (degrees) 


The aerodynamic forces and moments are simple linear rough 


estimates and do not take into account non-linearities or 


interaction of the aircraft with the store [Refs. 2, 6, 


Ze 


and 7]. As the design progresses past the preliminary stage a 


more detailed airload analysis must be performed. 


IO 
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V. COMPUTER PROGRAM 


A computer program to conduct an analysis of the loads at 
the aircraft-store interface in accordance with Mil-A-8591H, 
was written by the author and installed on the VAX/VMS 
Computer System, Advanced Computational Laboratory of 
Aeronautics and Astronautics Department, Naval Postgraduate 
School. The program is designed to be utilized by aero 
engineering graduate students in their aircraft design 
courses. The input can be interactive or batch. The 
interactive input is useful when running the program for the 
first time. The batch input is the most convenient when 
several runs are required or when the student requires a copy 
of their data input. 

The program flow chart is depicted in Figure 6. The 
program consists of an input subroutine, an output subroutine, 
an inertial load subroutine, an airload subroutine, a 
Subroutine to combine inertial loads and airloads, four 
Subroutines to calculate the loads at the respective 
attachment points, and a subroutine to resolve the store 
coordinate system with the aircraft coordinate system (for 
stores mounted at an angle to the aircraft coordinate system). 

Once the data is input to the program, the airloads (if 
applicable) on the store CG are calculated, and the inertial 


loads are calculated. The loads are then combined and the 


Ze 








Program Flowchart 


Figure 6. 
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appropriate coordinate system is selected. The loading on the 
store is analyzed. The critical attachment component is 
determined and the appropriate set of equations are utilized 
to determine the worst case load on the appropriate store 
attachment component. The resultS are sent to the computer 
screen and an output file. The input and output data are 
available to the user until the program is run again. The 
files contain the data of the most recent run. If hard copies 
are desired they must be made before rerunning the program. 
Detailed instruction on using the program is provided in the 


Mil-A-8591H Users Manual, Appendix A of this text. 
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VI. RESULTS AND RECOMMENDATIONS 


The program was tested against the sample problem in Mil- 
A-8591H, Appendix D [Ref. 1], of the specification. The sample 
consisted of a generic 2000 1b bomb analyzed with the general 
store inertial loading of Mil-A-8591H, Appendix A [Ref. 1]. 
The sample problem analyzed two inertial loading cases for the 
forward sway braces, two inertial loading cases for the aft 
Sway braces, two inertial loading cases for the forward 
hook/lug, and two inertial loading cases for the aft hook/lug. 
The results of the computer program were identical to the 
results provided in the sample problem. The only difference is 
the specification sample problem used five significant figures 
while the computer program uses eight significant figures. 

There was no airload analysis in the specification sample 
problem. The airload analysis in MIL-A-8591H is the same as 
the airload analysis in a current loads program [Ref. 2]. The 
inertial loads and airloads were combined and the resultant 
forces calculated. The effects of the airloads changed the 
resultant as expected. The program worked for all sample cases 
tested. 

The program was successfully incorporated into the Naval 
Postgraduate School Aero Department VAX computer. The results 
show the data to be identical to the sample data generated in 


Mil-A-8591H, Appendix D [Ref. 1). The computer code 
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accomplishes the objective of this thesis, to provide an 
analysis of the aircraft-store interface reactions in 
accordance with Mil-A-8591H. The program will be useful as a 
design tool in aircraft, helicopter, and weapons design 
Classes. Further testing and analysis of the program and 
methods incorporated in the program are required before this 
program is used in actual aircraft-store integration. This can 
be conducted, while the program is being used in aircraft, 
helicopter, and missile design courses at the Naval 
Postgraduate School. 

Recommendations include comparing the results of the 
MIL-A-8591H loads program with a finite element analysis of a 
specific store. The results indicate the validity of the 
assumptions made in deriving the equations for the analysis. 

The computer code can be made more efficient by creating 
a data base with all the store inertial envelopes of 
Mae A-8591H [Ref. 1], to calculate loads for the entire 
inertial envelope. From this, worst case load conditions for 
each sway brace and hook/lug can be identified and printed 
out. This will eliminate input errors and the possibility of 


the user not choosing the critical envelope loading condition. 


ie, 
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I. INTRODUCTION 


This user’s manual describes how to apply a computer 
program (MIL8591H) to calculate reaction loads on aircraft- 
Carried stores as recommended by military specification. The 
general deSign guide for aircraft store integration is 
Military Specification MIL-A-8591H, "Airborne Stores, 
Suspension Equipment And Aircraft-Store Interface (Carriage 
Phase); General Design Criteria For" [Ret. Hj. -This 
Specification outlines a method to calculate reaction loads 
for an aircraft suspension system equipment / store interface. 
A store is defined as any equipment that can be temporarily 
mounted to an aircraft wing or bomb bay pylon. Some common 
stores include weapons, auxiliary fuel tanks, electronic 
Surveillance pods, electronic countermeasure pods, and gun 
pods. Weapons make up the majority of store types. Weapon 
Suspension equipment is the interface of the store to the 
aircraft. Most stores are mounted to the aircraft by 
Suspension equipment consisting of two hook/lug assemblies, 
and further secured by two sets of sway brace pads. A standard 
configuration is illustrated in Figure 7. The Aircraft Stores 
Interface Manuals (ASIM) Volumes 1-3 [Refs. 3, 4, and 5] 
provide detailed descriptions and engineering drawings of 


stores (Vol. 1), Suspension equipment (Vol. 2), and aircraft 


Sil 


(Vol. 3). These references are currently utilized by the 


Department of Defense and NATO COUMEmaece 


SZ 


II. DESCRIPTION 


The purpose of the MIL8591H program is to calculate the 
worst case load for each mounting interface (i.e., forward 
hook/lug, aft hook/lug, the four sway braces forward left, 
mouwanscd right, aft left,aand aft right.) The program applies 
inertial and airloads to produce a net reaction at each 
attachment point. To successfully run the program three sets 
of information must be provided. The three major inputs are 
store/suspension equipment physical characteristics, inertial 
load factors and airload factors. A separate detailed section 
will be devoted to illustrate how to obtain this information. 
Figure 8 presented a flowchart that summarizes how the 


computer program works. 
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III. PHYSICAL CHARACTERISTICS 


The major physical characteristics are the stores weight, 
mass moments of inertia, the distance of the store center of 
gravity to the hook/lugs and the distance of the center of 
gravity to the sway braces. For most current weapons and 
Suspension equipment this information can be found in The 
Aircraft Stores Interface Manual (ASIM) Volumes 1-3; [Refs. 3, 
4a nd goa: 

Weight of the store, which will be abbreviated W in the 
computer program, and the Mass Moments of Inertia, which will 
be abbreviated IX, IY, and IZ, can be obtained directly from 
Volume 1 of the ASIM Manual, Store Characteristics [Ref. 3]. 
The weight must be entered in units of pounds and the mass 
moments of inertia must be entered in units of slug-ft%*2. Note 
that the program automatically converts the mass moments of 
inertia to units of lb-in-sec’* during calculation to stay 
compatible with English (U.S.) units in the equations 
developed by the specification. 

The distance of the center of gravity (CG) to the hook/lug 
and the distance of the center of gravity to the Sway braces 
must be calculated. The first step is to directly read the 
distances from the CG to the hook/lugs using Volume 2 of the 
ASIM Manual, Store Characteristics [Ref. 4]. All distances are 


to be read in inches and as positive. The computer program 
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uses only positive distances to calculate the loads. There are 
seven distances with respect to CG and hook/lug location that 
are of concern. The distances are referred to as inputs L, Ll, 
ime 5o,. b4, 5, andes (see Figures 9 and 10). The first 
input, the hook/lug spacing (L), is the distance from the 
forward store hook to the aft store hook. Inputs Ll and L2 are 
the longitudinal distances (x-direction) from the store CG to 
the forward hook/lug and to the aft hook/lug respectively. 
Inputs L3 and L4 are the lateral distance (y-direction) from 
the store CG to the forward and aft hook/lug respectively. The 
inputs L5 and L6 are the vertical (z-direction) distances from 
the store CG to the forward and aft hook/lug respectively. For 
most stores, L3 uSually equals L4, and LS equals L6 but 
caution should be exercised when dealing with non symmetric 
stores or stores mounted at an angle to the vertical. Figures 
9 and 10 illustrate the store physical characteristics with 
respect to the hook/lug. 

The distances of the CG to the sway braces are calculated 
using Volume 3 of the ASIM Manuals, Suspension Equipment 
[Ref. 5]. Detailed drawings of commonly used suspension 
equipment are presented. The problem is to relate the physical 
characteristics of the store to the physical characteristics 
of the suspension equipment. A point common to both the store 
and the suspension equipment is the hook/lug. By transforming 
the relation of the store CG with respect to the store lug to 


the suspension equipment drawing, the CG can be related to the 
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suspension equipment hook and the distances from the CG to the 
Sway braces can be obtained. 

There are seven inputs with respect to CG and sway brace 
location that are required. The inputs are referred to as S, 
S1, S2, S3, S4, S5, and S6 (see Figures 9 and 10). The first 
input, the sway brace spacing (S), is the distance from the 
forward sway brace to the aft sway brace. The inputs Sl and S2 
are the longitudinal distances (x-direction) from the CG to 
the forward sway brace and to the aft sway brace respectively. 
The inputs S3 and S4 are the lateral distances (y-direction) 
from the CG to the left and right sway brace respectively. The 
inputs S5 and S6 are the vertical (z-direction) distances from 
the CG to the forward and aft Sway brace respectively. For 
most stores, S3 uSually equals S4 and S5 equals S6, but 
caution should be exercised when dealing with non symmetric 
stores or stores mounted at an angle to the vertical. Figures 
9 and 10 illustrate the store physical characteristics with 
respect to the sway braces. 

The only other physical characteristics required to run 
the load analysis are the angle of the sway brace pads with 
respect to the vertical, designated (EPS) epsilon; and the 
angle of the mounted store with respect to the vertical, 


designated (DEL1) delta 1. Figure 10 illustrates both angles. 
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IV. INERTIAL LOADS 


The inertial loads required for the analysis consist of 
Six accelerations to the store’s center of gravity. They are 
the longitudinal acceleration (Nx), lateral acceleration (Ny), 
vertical acceleration (Nz), roll angular acceleration (PHIDD), 
pitch angular acceleration CHE DE) and yaw angular 
acceleration (PSIDD). Note that PHIDD, THEDD,and PSIDD are 
snort for phi double dot, theta double dot and psi double dot 
and correspond to the NASA notation for the angular 
acceleration of an air vehicle. 

The coordinate system and sign convention (right hand 
coordinate system) utilized in the program is illustrated in 
Baogure il. The sign convention of the accelerations is an 
important consideration in the analysis of the problem. All 
the signs of angles and rates must be consistent with the 
positive coordinate system indicated, Figure 11. An incorrect 
Sign will invalidate the entire analysis. The longitudinal 
axis (x-axis) 1S poSitive aft and roll acceleration is 
clockwise looking aft. The lateral axis (y-axis) is positive 
right looking down and the pitch acceleration is positive 
pitch nose up. The vertical axis (z-axis) is positive up and 
positive yaw acceleration iS positive nose left looking 


forward. 


ony 


General inertial envelopes are provided for tactical 
carrier based aircraft in Figures 12 and 13, and for 
helicopters in Figure 14. These inertial envelopes are to be 
utilized when no individual aircraft carriage is specified or 
when a broad spectrum of carriage aircraft are being 
considered. When more specific aircraft performance data is 
available the inertial envelopes can be adjusted accordingly 
with the consent of the acquiring activity. Mil-A-8591H, 
Appendix B [Ref. 1], can be utilized to give more specific 
envelopes. Since the analysis being conducted is for 
preliminary designs, the general envelopes will be utilized. 

The inertial envelope for the tactical carrier aircraft 
Figure 12, contains envelopes for arrested landings, catapult 
take offs, and in-flight wing-mounted stores. A second in- 
flight inertial envelope is presented for wing-tip mounted 
stores. The inertial envelope for the tactical carrier 
aircraft Figure 13, contains envelopes for arrested landings, 
catapult take offs, and in-flight fuselage mounted stores. 
Notice that the wing mounted stores have the most Severe 
inertial loading. The inertial envelopes in Figures 12, 13, 
and 14, are actually three dimensional envelopes, Since there 
is aeplus/minus g factor listed with the angular 
accelerations. 

To utilize the computer program, a worst case inertial 
loading must be determined. The worst case inertial loadings 


will be found at the "corners" of the inertial envelopes. Mil- 
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A-8591H, Appendix D [Ref. 1], describes a procedure to 
determine the worst case load for each suspension equipment 
component (i.e. hook/lug or Sway brace). By assuming that the 
hook/lug can only take a tensile load and the sway braces can 
only take a compressive load the direction of the worst case 
inertial loading is presented in Figure 15. For example, the 
worst case inertial condition for the forward left sway brace 
is when there is a negative x acceleration, a negative y 
acceleration, a positive z acceleration, a negative x angular 
acceleration (roll), a positive y angular acceleration 
(pitch), and positive z angular acceleration (yaw). Note that 
the sign conventions of Figure 11 must be adhered to for a 
valid computer analysis. Once the direction of the worst case 
inertial loads are determined, the magnitudes are obtained 
from the envelopes depicted in Figures 12 and 13. The computer 
program automatically indicates to the user the component that 
will experience the worst load when an inertial loading 
condition is applied. 

To utilize the analysis for a helicopter, the directions 
of the acceleration are determined by Figure 15 while the 
magnitudes of the acceleration are determined by Figure 14. 
The inertial envelope produced in Figure 14 is a generic 
envelope, calculated using generic data from the equations of 
Mil-H-8591H, Appendix C [{Ref. 1]. The inertial load envelopes 
of Figure 14 are only for preliminary design. Specification 


compliance with this specific helicopter inertial load 
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envelope is not mandated but requires the approval of the 
aCqulringeactaswaty . 

The reaction loads at the store interface are equal and 
opposite for the store and the suspension equipment, therefore 
the analysis can be used for preliminary store or aircraft 


design. 
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V. AIRLOADS 


Preloads are probably —~the most difficult forces to 
determine in the analysis. There are three options available 
to determine the airloads on the store. The first option is to 
Simply ignore the airloads. This is not a bad approximation 
for arrested landings and catapult takeoffs, but it is not 
realistic for in-flight data for high performance aircraft. 
The second option is to directly input the airloads (if 
available) from experimental data or wind tunnel modeling. 
Although a best estimate, this data is usually not available 
in the preliminary design stage. The third method is to 
Calculate the airloads from the data available with respect to 
maximum angles of attack, and maximum angle-of-sideslip. 

Mil-A-8591H, Appendix A [Ref. 1], provides equations to 
Calculate the maximum angle-of-attack and the maximum angle- 
of-sideslip for very specific flight maneuvers (i.e., 
pushovers, and pullups) at specific corners of the inertial 
flight envelope. The maximum angle-of-attack for the store 
(ANGA) and the maximum angle-of-sSideslip (ANGS) are functions 
of dynamic pressure and flight maneuver and can be calculated 
from Figure 16 for wing mounted stores and Figure 17 for 
fuselage mounted stores. The slope of the lift coefficient, 
Slope of the pitching moment coefficient, slope of the side 


force coefficient, and coefficient of drag must also be 
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available to perform the airload analysis. Frequently, these 
data are not available for a preliminary design. The method 
presently applied [Ref. 2] models the store as a simple 
cylinder (most Stores are) and usesS generic coefficients to 
make an approximate of the airloads [Refs. 6, 7]. The computer 
program uses this option as a default (model the store as a 
Simple cylinder) capability. See the sample problem for 
implementation of the airload default capability. 
Computation of airloads for helicopters is not addressed 
in the specification. The recommended method of helicopter 
airload analysis (i.e., wind tunnel analysis or rotor flight 
Simulation programs) is beyond the preliminary design phase. 
Considering the lower airspeed that helicopters fly, the 
initial estimate of airloads with respect to inertial loads is 
not aS critical. Unless the airloads of a store are available, 
the next best option is to estimate the airloads during the 


preliminary design stage. 
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VI. ASSUMPTIONS AND LIMITATIONS 


Due to the fact that the store is ae statically 
indeterminate structure with unknown flexibility 
characteristics, the analysis to calculate the reaction on the 
store suspension equipment is a relatively complex problem. By 
making several assumptions, the problem can be simplified so 
a general analysis can be performed for stores and aircraft 
that are still in the preliminary design stage. 

The analysis covered in this program assumes that the 
store is a rigid body. Other assumptions are: the loads at the 
hook/lug are tensile only, and the loads at the sway brace 
pads are compressive only; the sway braces are usually at a 
roll angle to the store and can exert a lateral and vertical 
force, but cannot exert a longitudinal force. Since designers 
are concerned with worst case loads the equations are derived 
using the external load conditions of Figure 15. These 
external loads provide the maximum reaction forces at each 
individual hook/lug and sway brace as indicated. The result of 
these assumptions simplifies the problem to a statically 
determinate problem. The analysis includes both inertial and 
airloads. Other contributions to load factors such as flutter, 
rotor down wash, modal coupling, are beyond the scope of this 


analysis and beyond the scope of the preliminary design stage, 
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but must be addressed during the detail design Stage and 


before the critical design review (Cie 
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VII. UTILIZING THE PROGRAM 


The MIL8591H fortran computer program is installed on the 
VAX/VMS Computer System, Advanced Computational Laboratory of 
Aeronautics & Astronautics Department, Naval Postgraduate 
school. The program consists basically of three parts; input 
file, output file, and program code. The input file is called 
eeorUl0.dat", the output file is called “for0ll.dat"”, and the 
program code is called "Mil8591H.for". 

A description and implementation of the files is described 


in the following paragraphs. 
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VIII. INPUT DATA 


The first step toward executing the program is to input 
the data. Data input can be either interactive or by a batch 
input file. The interactive data input is self explanatory but 
time consuming when multiple runs are required. A batch file 
is efficient to run several design cases. A hard copy of the 
batch input file could also be printed for later study. A 
first run with interactive is recommended, but once the user 
becomes familiar with the program, batch inputs are more 
convenient and efficient. 

To demonstrate the execution of the program, a batch file 
will be created. For the program to work, the batch input file 
must be named, for010.dat. The batch input file will already 


exist in the computer and will look like the following: 


STORE CG POSITION WITH RESPECT TO SWAY BRACE PADS. 

Daas CG POSITION Wilh RESPECT, 1O HOor, sUGs- 

ae WEIGHT (W) LBS 

1000 

STORE MOMENT OF INERTIA X-AXIS (IX) SLUG-FT**2 

er MOMENT OF INERTIA Y-AXIS W1Y) poLUG-E 1]. Z 

oe MOMENT OF INERTIA Za xiS)(i2) )cuUG=-n 1-2 

ae BETWEEN FOR AND AFT SWAY BRACE PADS (S) IN 
eee BETWEEN CG AND FOR SWAY BRACE PADS X-DIRECTION (S1) 
ae 
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DISTANCE BETWEEN CG AND AFT SWAY BRACE PADS X-DIRECTION (S2) 
en. 

9 

DISTANCE BETWEEN CG AND NEAR SWAY BRACE PADS Y-DIRECTION (S3) 
IN. 

Bie 2 

DISTANCE BETWEEN CG AND FAR SWAY BRACE PADS Y-DIRECTION (S4) 
IN. 

Sey 

DISTANCE BETWEEN CG AND RIGHT SWAY BRACE PADS Z-DIRECTION (S5) 
oN 

oe 0 

DISTANCE BETWEEN CG AND LEFT SWAY BRACE PADS Z-DIRECTION (S6) 
ILE 

oa0 

DISTANCE BETWEEN FOR AND AFT LUGS X-DIRECTION (L) IN. 

14 

DISTANCE BETWEEN CG AND FWD LUG X-DIRECTION (L1) IN. 

5 

DISTANCE BETWEEN CG AND AFT LUG X-DIRECTION (L2) IN. 

5 

DISTANCE BETWEEN CG AND FWD LUG Y-DIRECTION (L3) IN. 

0 

DISTANCE BETWEEN CG AND AFT LUG Y-DIRECTION (L4) IN. 

0 

DISTANCE BETWEEN CG AND FWD LUG Z-DIRECTION (L5) IN. 

ia) 

DISTANCE BETWEEN CG AND AFT LUG Z-DIRECTION (L6) IN. 


erat LOAD FACTOR X-DIRECTION (NX). 
aa LOAD FACTOR Y-DIRECTION (NY). 
teria. LOAD FACTOR Z-DIRECTION (NZ). 
EXTERNAL ANGULAR ACCELERATION X-DIRECTION (PHIDD) RAD/SEC%2. 


EXTERNAL ANGULAR ACCELERATION Y-DIRECTION (THEDD) RAD/SEC%2. 
-4 

EXTERNAL ANGULAR ACCELERATION Z-DIRECTION (PSIDD) RAD/SEC%2. 
=Z 

ANGLE BETWEEN RADIUS OF CURVATURE OF STORE AND Z-AXIS (EPS) 
DEG. 

20 

THE ROLL ANGLE BETWEEN THE AIRCRAFT AND THE STORE (DEL1) DEG. 
0 

HOW AirloadS ARE TO BE ANALYZED. ENTER OPTION 1,2,OR 3. 

3 

Soele THE ATRLOADS ARE NEGLECTED OPTION 1 DO NOT ENTER ANY 
FURTHER DATA 
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Cx IF THE AIRLOADS ARE TO BE ANALYZED (I.E. OPTION 2 35pe 
selected) 

Cx OPTION 2 THE FOLLOWING 6 LINES MUST BE ENTERED 

C* OPTION 3 THE LAST 12 LINES MUST BE ENTERED 


C* INPUT AIRLOADS DIRECTLY FROM WIND TUNNEL DATA OPTION 2 


LONGITUDINAL AIRLOAD (PXA) LB X-DIRECTION POSITIVE AFT 
ease AIRLOAD (PYA) LB Y-DIRECTION POSITIVE Lier 
VERTION. AIRLOAD (PZA) LB Z-DIRECTION POSITIVE UP 

aT MOMENT- ROLL (AMXA) IN-LB POSITIVE CCW LOOKING FORWARD 
eves MOMENT— PITCH (AMYA) /IN-LE POS Ret Eaves raor 

va MOMENT- YAW (AMZA) IN-LB POSITIVE NOSE LEFT 


Gx 

C*, INPUT DATA TO CALCULATE ATRIEGADS. CellOne. 
Cx 

AIRCRAFT VELOCITY (VEL) - FT/SEC 

600 

ANGLE OF ATTACK (ANGA) - DEG 

76 

ANGLE OF SIDE-SLIP (ANGS) - DEG 

6 

AIR DENSITY (DEN) - SLUG/CUBIC FT 

-000Z7s 

REFERENCE RADIUS (R)- IN 

8 

C* WRITE STORE AERODYNAMIC CHARACTERISTICS (FOR UNKNOWN 
VALUES, 


C* ENTER A "0" AND DEFAULT VALUES BASED ON A CYLINDRICAL 
C* SHAPE WILL BE ENTERED 


REFERENCE VAREA. (SR) — FI so 

— EENGTH (be) - iN 

see LIFT: COBFPFICIENT SLOCE (Caer eo ohh rG 

on PITEGH MOMENT COEPEICIBNT SLOPE (C2) = fh, DEG 
— SIDE=FORCE COBFPICIENISSEGEE (G6) —" PER DEG 
STORE YAW MOMENT COEFFICIENT SLOPE (C7) - PER DEG 
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mloOkmE DRAG CORFFICIENT (CD) 
0 


The general format of the input file must remain as it is, 
only the input of numerical values can be changed. If comment 
lines are deleted or altered, it may not allow the program to 
Pain. 

tae. Hatch anput file inputs “STORE CG POSITION WITH 
meoemel TO SWAY BRACE PADS," "STORE POSITION WITH RESPECT TO 
HOOK/LUGS," and "HOW ARE THE AIRLOADS TO BE ANALYZED" need 
explanation. All other inputs are self explanatory. 

ine input "STORE CG POSITION WITH RESPECT TO SWAY BRACE 
PADS" requires a number from 1 to 7 to be entered. The 
following cases represent the options available. 


ie The store CG is between the sway braces longitudinally 
and laterally. 


2 The store CG is between the sway brace pads 
longitudinally but outside the sway brace pads 
laterally. 

Bis The store CG is forward of the front sway brace pad 
longitudinally but between the sway brace pads 
laterally. 

4. The store CG is forward of the front sway brace pad 
longitudinally but outside the sway brace pads 
eee re celle 

as The store CG is between the sway brace pads 
longitudinally but outside the sway brace pads 
laterally. 

og The store CG is aft of the rear sway brace pad 
longitudinally but between the sway brace pads 
laterally. 
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as The store CG is aft of the rear sway brace pad 

longitudinally but outside the Sway brace _ pads 
laterally. 

The large majority of stores fits into category 1 with the 

store CG located between the sway braces longitudinally and 

laterally. 
The input "STORE CG POSITION WITH RESPECT TO HOOK/LUGS"” 
requires a number from 1 to 3 to be entered. The following 


cases represent the options available. 


1 The store CG is between the hook/lugs. 
es The store CG is in front of the forward hook/lug. 


ce The store CG is aft of the rear hook/lug. 


The large majority of stores fits into category 1 with the 
store CG located between the hook/lugs. 

The input "HOW ARE THE AIRLOADS TO BE ANALYZED" requires 
a number from 1 to 3 to be entered. The following cases 


represent the options available. 


jigs Neglect the airloads (i.e. arrested landing, catapult 
takeoff) 

Ze Input airloads directly from wind tunnel data. 

ce Calculate airloads. 


Option 1 neglects the airloads and assigns values of zero 
for the longitudinal, lateral, and vertical airloads. The 
aerodynamic roll pitch and yaw moment are also assigned a 


value of zero. Option 1 would be used when an analysis of only 
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inertial loads are required as in arrested landings or 
catapult takeoffs. If Option 1 is entered the program stops 
reading the input file at that data point and proceeds with 
calculations. The results display the zero airload condition. 

Option 2 allows the user to directly input airloads and 
aerodynamic moments obtained by other methods such as wind 
tunnel data or data from Similar stores. If Option 2 is 
entered the six data points PXA (longitudinal airload), PYA 
(lateral airload), PZA (vertical airload), AMXA (aerodynamic 
roll moment), AMYA (aerodynamic pitch moment), and AMZA 
(aerodynamic yaw moment) will be read into the program and 
combined with the inertial loads to calculate the net reaction 
forces at the critical sway brace or hook/lug. The output will 
display the aerodynamic loads and moments. 

Option 3 will calculate the aerodynamic loads and moments. 
There are 12 inputs required to calculate the loads and 
moments. The inputs VEL (aircraft velocity), ANGA (store angle 
of attack), ANGS (store angle of sideslip), DEN (air density), 
and R(store reference radius) must be entered, since no 
default values are available for these five inputs. The 
maximum angle-of-attack for the store (ANGA) and the maximum 
angle-of-sideslip (ANGS) are functions of dynamic pressure and 
flight maneuver and can be calculated from Figure 16 for wing 
mounted stores and Figure 17 for fuselage mounted stores. The 
final six inputs; SR (reference area), LR (reference length), 


Cl (store lift coefficient slope), C2 (store pitch coefficient 


oH 


slope), C6 (store side-force coefficient slope), C7 (store yaw 
moment coefficient slope), can be entered directly if the 
values are known or a zero can be entered and default values 
based on a cylindrical shape will be utilized in the 
calculation of the aerodynamic forces and moments. The airload 
analysis based on the cylindrical shape [Refs. 6 and 7] is 
utilized in a current MIL-A-8591 loads analysis [Ref. 2]. The 
output will display the aerodynamic loads and moments used in 


calculation of the reaction forces on the store. 
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IX. RUNNING THE PROGRAM 


Assuming the user logged on to the computer and edited the 

input file for010.dat with the desired inputs simply type; 
run mMieleooih 
The program will prompt the following; 
ENTER 1 FOR INTERACTIVE 2 FOR BATCH I0 

Enter the desired number and the program will continue. If you 
run interactive you will enter the same data, one line at a 
time, as in the batch input file. Once the data is read, the 
program automatically runs and output data is displayed on the 
screen and saved in output data file, for0Oll.dat. To obtain a 
hard copy of the input file and/or the output file print the 
appropriate file. For example on the NPGS AERO DEPARTMENT VAX: 

To obtain a hard copy of the input file type; 

print for010.dat/que=1n03 
To obtain a hard copy of the output file type; 
print for011.dat/que=1n03 

The input file will look like the example in the input 
section of this manual and the output file will look like the 
example in the output section of this manual. If hard copies 
are desired they must be made prior to rerunning the program. 
The input and output file contain only the latest data i.e., 
the input and output files are updated with each run. To run 


the program again, edit the input file, for010.dat, with the 


Se 


desired changes in conditions and Cype run mil859ih. ae 


input and output file will be created with the latest data. 
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X. OUTPUT FILE 


The output data will indicate which component (lug or sway 
brace) experiences the worst case load. The magnitude of the 
loads for the indicated directions will be presented. Remember 
the hook/lugs can only react in tension and the sway braces 
can only react in compression. The Sway braces cannot react a 
load in the longitudinal or x-direction. The inertial loads 
and airloads applied in the analysis are presented. The major 
physical characteristics, weight and mass moments of inertia 
are indicated, and finally the data utilized to compute the 
airloads, if applicable, are presented. Since physical 
characteristics usually remain constant for a store, a record 
of analysis for a particular store would consist of one input 
file attached to the multiple output files. 

The actual output file for the data represented in the 


input section of this manual will look like the following; 


WORST CASE LOAD IS AT THE AFT HOOK 


AFT HOOK/LUG LOADS LBS 


VERTICAL = 18977.46 
LATERAL = 659.915 90 
HONGITUDINAL = 1290.393 
TOTAL SHEAR = 1449.379 


*XINERTIAL LOADS APPLIED** 
LOAD FACTOR G"S ANGULAR ACCELERATION RAD/SEC**2 


LONGITUDINAL = -1.500000 ROLL = 0.0000000E+00 
LATERAL =. =5. 000000 PEten = -4.000000 
VERTICAL = -12.00000 YAW = 2 000000 
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**SAERODYNAMIC LOADS APE ir D a> 


AIRLOAD LB AERODYNAMIC MOMENT IN-LB 
LONGITUDINAL = 20926060 ROLL = 0.0000000E+00 
LATERAL = 134.1484 PRotch = P21 49326 
VERTICAL = 699224 YAW = Se AUD a a1 
STORE WEIGHT LBS = 2000.000 
MASS MOMENT OF INERTIA X-AXIS LB-IN-SEC*2 = 57.60000 
MASS MOMENT OF INERTIA Y-AXIS LB-IN-SEC*%2 = 1272.000 
MASS MOMENT OF INERTIA Z-AXIS LB-IN-SEC°2 = 127223000 


HOOK/LUG SPACING IN. = 14.00000 
THE STORE CG IS 7.00000 IN PROM THE Fuh Ook fuc 


LONGITUDINAL SWAY BRACE SPACING IN = 18.00000 
THE STORE CG IS 9.00000 IN FROM THE FWD SWAY BRACE 


AIRLOADS WERE CALCULATED CASE 3 


ATVRERAY | VigiOC Diese / SiC = 600.0000 
STORE ANGLE OF ATTACK DEG = 76.00000 
STORE ANGLE OF SIDESLALP DEG & ©.0000uy 
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XI. SAMPLE PROBLEM 


The purpose of the sample problem is to demonstrate 
application of the computer program for a generic store and 
generic suspension equipment (i.e., bomb and bomb rack). The 
sample configuration is depicted in Figure 18. The necessary 
dimensions to run the program are indicated. The computer 
program can be utilized either interactively or by a batch 
mmelei. Lf multiple iterations are to be analyzed, it is 
recommended to use the batch input. ReSultS are provided on 
the user’s screen and the results are recorded in an output 
file so a hard copy can be produced. A summary of changes to 
the input file with a summary of the resulting changes in the 
output file will illustrate execution of the problem. The 
details of operating the computer are in the section of this 
manual on the running the program. 

The analysis is accomplished uSing three steps. The first 
step is to enter the physical characteristics into the input 
rile. 

The second step is to input an inertial load from the 
appropriate inertial envelope Figures 12, 13, or 14, and 
critical load Figure 15. Rerun the program with different 
inertial load conditions until you are satisfied you have the 
most critical case. The third step is to apply the appropriate 


mmerhoads to the most critical *inflight* inertial load 


oe 


condition. By comparing the catapult, arrested, and in-flweme 
(with airloads) cases, the worst case load can be identified. 

The physical characteristics of the store and suspension 
equipment are illustrated in Figures 9 and 18. Edit the input 
file FOROQ10.DAT and replace the appropriate data. The 
respective input file is depicted in the input section of this 
manual. Enter the inertial load condition and initially select 
"1" (neglect airloads) under the "HOW DO YOU WANT THE AIRLOADS 
ANALYZED" input. 

This allows the user to find the worst case inertial load. 
After the worst case inertial loads are calculated, add 
airloads to the in-flight case by rerunning the program with 
a "2" or "3" entered under the "HOW DO YOU WANT THE AIRLOADS 
ANALYZED" input. The appropriate aerodynamic data must also be 


entered. 


A. FRONT SWAY BRACE 

In this example we analyze the left front sway brace for 
a wing mounted store on a carrier based aircraft. From Figure 
15 the critical store inertial load for the left front sway 
brace 1S —-NX, “NY, +NZ,7 =PHIDD, 37 2eebl ane +PSIDD. De 
envelope in Figure 12 is utilized for wing mounted stores on 
a carrier based aircraft to find the magnitudes at the worst 
case. In this example there are several possibilities for a 


worst case (corners of the envelope) loading. 
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Remember to keep sign conventions consistent as indicated. 


NX NY NZ ) W ENVELOPE 
@ASE 1: -9 -5 +2 -100 +25 +6 ARRESTED LANDING 
CASE 2: -3 -5 +4 -100 +25 +6 ARRESTED LANDING 
Saor 3: isos Vado, he 0 +4 +2 FLIGHT 


By running the program with these store inertial loads the 


results for the forward left sway brace are; 


SZ) (Lbs) Sy (lbs) Seoutlps) ENVELOPE 
CASE 1: LOZ 315). 5 9) Zs ARRESTED LANDING 
Ao 2: 9636 8y53) 7h 10254 ARRESTED LANDING 
GASE 3: 4150 10 4416 EEIGHT 


Now CASE 3 must be analyzed with airloads added. Since the 
component is the left forward sway brace, the critical 
airloads will be in the positive vertical (z-direction), 
negative lateral (y-direction), positive store pitch moment 
(nose up), and positive store yaw moment (nose left). If 
option 2 for analyzing the airloads is selected the data is 
entered with the sign and magnitude under the appropriate 
entry. 

If option 3 is selected, Figure 16 is utilized to find the 


largest angle-of-attack and largest angle-of-sideslip to 


Sa, 


produce the desired airload signs. In this example a positive 
angle-of-attack and positive angle-of-sideslip will produce 
the critical airloads. To find the magnitude of these angles, 
assume that the aircraft is travelling at 600 ft/sec at an air 
density of .00278 slugs/ft’?. Using the equation for point 4 
(the in-flight inertial point of CASE 3), the maximum angle- 
of-attack is 0 degrees and maximum angle-of-sideslip is +6 
degrees. Enter the inputs into the input file (in this example 
enter default values "0" for the coefficient slopes) and rerun 
the program. The result for CASE 3 with airloads is 5779 lbs. 
The worst case load for the left front Sway brace is CASE 1 
ARRESTED LANDING. 

The worst case load for the right front sway brace is the 
Same as the left front sway brace, since the magnitudes of the 
angular accelerations are equal in the positive and negative 
direction. The store can be mounted on the left or right wing 
cauSing the in-flight envelope, Figures 12 and 13, to be 


symmetric. 


B. AFT SWAY BRACE 

In this example we analyze the left aft sway brace for a 
wing mounted store on a carrier based aircraft. From Figure 15 
the critical store inertial load for the left front sway brace 
is +NX, -NY, +NZ, -PHIDD, -THEDD, -PSIDD. The envelope in 
Figure 12 is utilized for wing mounted stores on a carrier 


based aircraft to find the magnitudes at the worst case. In 
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this example there are several possibilities for a worst case 
(corners of the envelope) loading. 


Remember to keep sign conventions consistent as indicated. 


NX NY NZ d 6 W ENVELOPE 
CASE 1: +9 -2.5 +1.5 -25 -15 -4 CATAPULT T/O 
CASE 2: +2 -5 +4 -100 -25 -6 ARRESTED LANDING 
GSE 3: 41.5 -1.5 +7 0 -4  -2 FLIGHT 


By running the program with these store inertial loads the 


results for the aft left sway brace are; 


Sy a) ene) sy(lbs) Stot(lbs) ENVELOPE 


ease 1: S502 Zoe5 627 3 CATAPULT T/O 
Ao 2: 9401 3421 10005 ARRESTED LANDING 
SA5E 3: 4150 15 20 4416 PLIGHT 


Now CASE 3 must be analyzed with airloads added. Since the 
component is the left aft sway brace, the critical airloads 
will be in the positive vertical (z-direction), negative 
lateral (y-direction), negative store pitch moment (nose 
down), and negative store yaw moment (nose right). If option 
2 for analyzing the airloads is selected, the data is entered 


with the sign and magnitude under the appropriate entry. 
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If option 3 is selected, Figure 16 is utilized to fim ere 
largest angle-of-attack and largest angle-of-sidesiae ae. 
produce the desired airload signs. In this example a negative 
angle-of-attack and negative angle-of-sideslip will ppoGmee 
the critical airloads. To find the magnitude of these angles, 
assume that the aircraft is travelling at 600 ft/sec ac Mammen 
density of .00278 slugs/ft’*. Using the equation for jegmmee 
(the in-flight inertial point of CASE 3) maximum ne@amage 
angle-of-attack is -45 degrees and maximum negative angle-of- 
sideslip is -6 degrees. Enter the inputs into the impa@e ae 
(in this example enter default values "0" for the coefficient 
Slopes) and rerun the program. The results for CASE 3 Malena 
airloads is 8294 lbs. The worst case load for the left front 
Sway brace is CASE 2 ARRESTED LANDING. 

The worst case load for the right aft sway brace deme 
Same as the left aft sway brace since the magnitudes of the 
angular accelerations are equal in the positive and negative 
direction. The store can be mounted on the left or right wing 
causing the inflight envelope, Figures 12 and 13, tome 


symmetric. 


C. FORWARD HOOK/LUG 

In this example we analyze the forward hook/lug for a wing 
mounted store on a carrier based aircraft. From Figure 15 the 
critical store inertial load for the forward hook/lug is +NxX, 


+NY, -NZ, +PHIDD, -THEDD, +PSIDD. The envelope am Pigume aes 
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utilized for wing mounted stores on a carrier-based aircraft 
to find the magnitudes at the worst case. In this example 
there are several possibilities for a worst case (corners of 
the envelope) loading. 


Remember to keep Sign conventions consistent as indicated. 


NX NY NZ d 8 Ww ENVELOPE 
CASE 1: +9 -2.5 -5 -25 -15 -4 CATAPULT T/O 
CASE 2: +2 -5 -12 -~100 -25 -6 ARRESTED LANDING 
CASE 3: +1.5 -5 -12 0 -4  -2 FLIGHT 


By running the program with these store inertial loads the 


results for the forward hook/lug are; 


Lz (lbs) Lx(lbs) Ly(lbs) Lshear(lbs) ENVELOPE 


Seoemets 11382 9000 409 9009 CATAPULT T/O 
CASE 2: 15864 2000 Obs Da Ms: ARRESTED LANDING 
CASE 3: 12414 1500 678 1646 SMES ee 


Now CASE 3 must be analyzed with airloads added. Since the 
component is the left aft sway brace the critical airloads 
will be in the negative vertical (z-direction), negative 
lateral (y-direction), negative store pitch moment (nose 


down), and negative store yaw moment (nose right). If option 
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2 for analyzing the airloads is selected, the data is entered 
with the sign and magnitude under the appropriate entry. 

If option 3 is selected, Figure 16 is utilized to find the 
largest angle-of-attack and largest angle-of-sideslip to 
produce the desired airload signs. In this example a negative 
angle-of-attack and negative angle-of-sideslip will produce 
the critical airloads. To find the magnitude of these angles, 
assume that the aircraft is travelling at 600 ft/sec at an air 
density of .00278 slugs/ft’*. Using the equation for point 2 
(the inflight inertial point of CASE 3) maximum negative 
angle-of-attack is 0 degrees and maximum angle-of-sideslip is 
+6 degrees. Enter the inputs into the input file (in this 
example enter default values "0" for the coefficient slopes) 
and rerun the program. The results for CASE 3 with airloads is 
14111 lbs. The worst case load for the forward hook/lug is 


CASE 2 ARRESTED LANDING. 


D. AFT HOOK/LUG 

In this example we analyze the aft hook/lug for a wing 
mounted store on a Carrier based aircraft. From Figure 15 the 
Critical store inertial load for the aft hook/lug is -NX, #+NY, 
-NZ, +PHIDD, +THEDD, +PSIDD. The envelope in Figure 12 is 
utilized for wing mounted stores on a Carrier based aircraft 
to find the magnitudes at the worst case. In this example 
there are several possibilities for a worst case (corners of 


the envelope) loading. 


64 


Remember to keep sign conventions consistent as indicated. 


NX NY NZ od 6 WV ENVELOPE 
CASE 1: -3 -5 -12 -100 -25 -6 ARRESTED LANDING 
CASE 2: -9 -5 -6 -100 -25  -6 ARRESTED LANDING 
CASE 3: -1.5 -5 -12 0 -4 -2 FLIGHT 


By running the program with these store inertial loads the 


results for the forward hook/lug are; 


Lz (lbs) Lx(lbs) Ly(lbs) Lshear(lbs) ENVELOPE 


Stee i: 16332 3000 956 3149 ARRESTED LANDING 
CASE 2: 16145 9000 958 9050 ARRESTED LANDING 
CASE 3: 12414 1500 678 1646 FEIGHT 


Now CASE 3 must be analyzed with airloads added. Since the 
component is the left aft Sway brace the critical airloads 
will be in the negative vertical (z-direction), negative 
lateral (y-direction), positive store pitch moment (nose up), 
and negative store yaw moment (nose right). If option 2 for 
analyzing the airloads is selected the data is entered with 
the sign and magnitude under the appropriate entry. 

If option 3 is selected, Figure 16 is utilized to find the 
largest angl-of-attack and largest angle-of-sideslip to 


produce the desired airload signs. In this example a positive 
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angle-of-attack and negative angle-of-sideslip will produce 
the critical airloads. To find the magnitude of these angles, 
assume that the aircraft is travelling at 600 ft/sec at an air 
density of .00278 slugs/ft*. Using the equation for pommeme 
(the inflight inertial point of CASE 3) maximum positive 
angle-of-attack is 76 degrees and maximum angle-of-sideslip is 
+6 degrees. Enter the inputs into the input file (in this 
example enter default values "0" for the coefficient slopes) 
and rerun the program. The results for CASE 3 with airloads is 
18977 lbs. The worst case load for the aft hook/lug is CASE 3 


FLIGHT WITH AIRLOADS. 
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Figure 7. Relationship of Typical Bomb-Rack Top-View 
Drawing From Section 1 to the Actual Rack and the 
Carrying Aircraft. 
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Set Alrioads = 0 





s A 
oF: 
¥ s 
i ii BH 


i 


F 





Program Flowchart 


Figure 8. 
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Figure 9. Store Physical Characteristics 
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STORE CG 





Y-Z PLANE LOOKING FORWARD 


DEL1=0 for the analysis of the hook/lug indicated 


EPS 









STORE CG 


SWAY BRACE 


Y-Z PLANE LOOKING FORWARD 
DEL1#0 for the analysis of the hook/lug indicated 


Figure 10. Store Physical Characteristics and Angles 
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Load-Factor Envelope for 
Wing Mounted Stores 





Figure 11. Coordinate System, Sign Convention, and a Typical 
Load Factor Envelope 


fis 





CATAPULTING 


= 225 
= + 15 RAD/SEC 


= + 4RADSEC 
= + 25 RADSEC 





Figure 12. Design Inertial Limit Load Factors for Wing or 
Sponson-Mounted Stores (Data Applies at the 
Center of Gravity) 
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Figure 13. Design Inertia Limit Load Factors for Fuselage- 
Mounted Stores (Data Applies at the Store Center 
of Gravity) 


qs 


FLIGHT WING OR SPONSON 
MOUNTED STORES 


HELICOPTERS 

Pe Eels 6.0-1.5 

6 «= + 3 RAD/SEC# 
¥ = + 2 RAD/SEC# 


n 





6 = pitch acceleration 


FLIGHT FUSELAGE 
MOUNTED STORES 


HELICOPTERS 

n, = +6.0-1.5 

6  =+0.5 RAD/SEC2 
# =+1RAD/SEC 





Figure 14. Inertial Load Factors for Helicopters 
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SEQre Inertial SB/store SB pad moor ue 
Loads 
Fwd Aft 
L L A 
= + 


el A Oe Ee A 
| 
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NOTE: (1) If the CG of the store is located laterally outside 
of the SB’s, a negative vertical load (P,, down) may be 
Critical and should be investigated. 


Figure 15. Direction of External Loads and Moments for 
Maximum Reaction Forces at the SB’s/Stores SB 
Pads and Hooks/Lugs 


qh, 


For all points, the stores shall be considered to be mounted at incidence angles of 
0 or -3 degress, whichever is more critical in each case, to be added to the values 
of a. given below: 


Points (1) and (2) (symmetric pulup): 
a, =O0b+ me degrees 
p sa degrees 


Points (3) and (4) (pushover): 


a =(0p. 2 degress 
b 4 

6 =+ SM depress 

_ 4 


Paint (5) (roling pushover) 
10 ». 1m 109.) degrees 


c= ae 


“oP degrees 
Point (6) (roling pullout) 
«On. Dane e dagrees 
p =+ HERD dopeee 


where 


INBOARD 





a. = Store angle of attack 


B, = store angie of sidesiip 


OY “eed 
q = > , dynamic pressure 


Figure 16. 


Store Angles of Attack and Sideslip at Specific 
Load Envelope Points for Wing or Sponson-Mounted 
Stores. 
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For all points, the stores shall be considered to be mounted at incidence angles of 
0 or -3 degress, whichever is more critical in each case, to be added to the values 
of a. given below: 


Points (1) and (2) (pullup): 
a. = 010+ degrees 
pA 


Points (3) and (4) (pushover: 
a “Ow. am degrees 
B. aw degrees 


qe 
a. 


OUTBOARD - LEFT OUTBOARD - RIGHT 


where 

a. = Store angie of attack 

B, = store angie of sideslip 
q = MY dynamic pressure 


DOWN 


Figure 17. Store Angles of Attack and Sideslip at Specific 
Load Envelope Points for Fuselage-Mounted Stores. 
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Vee STORE WEIGHT: -1900LBS 
i IX: 4.8 SLUG-FT2 
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LOOKING FORWARD 
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Figure 18. Sample Problem Store 


78 


DELI 


DEN 


Deas, 


L5 


L6 


LR 


XII. GLOSSARY 
Aerodynamic roll moment (in-1bs) 
Aerodynamic pitch moment (in-1lbs) 
Aerodynamic yaw moment (in-1bs) 
Store angle of attack of store (deg), @, 
Store angle of side-slip (deg), 6, 
Rotation around the X-axis for stores mounted at a 
roll angle W.R.T. the aircraft axis system. (deg) 
Pore density (Slug/ tes) 
Angle between radius of curvature of store and 
Z-axis (eps) deg. 
Mass moment of inertia X-axis slug-ft**2 
Mass moment of inertia Y-axis slug-ft**2 
Mass moment of inertia Z-axis slug-ft**2 
Distance between lugs (in) 
Distance store CG to fwd lug X-direction (in) 
Distance store CG to aft lug X-direction (in) 
Distance store CG to fwd lug Y-direction (in) 
Distance store CG to aft lug Y-direction (in). Note 
L3 usually equal to L4 
Distance store CG to fwd lug Z-direction (in) 
Distance store CG to aft lug Z-direction (in). Note 
L5 usually equal to L6 


Reference length (in) 
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LTOTAFT 


LTOTFEWD 


LXAF'T 


LXFWD 


ran 


LYFWD 


LZAFT 


LZFWD 


MX 


MXC 


MY 


MyC 


MZ 


MZC 


NX 


NY 


NZ 


PHIDD 


ESIDD 


Ex 


PXA 


Aft lug shear force X-Y plane (lbs) 

Forward lug shear force X-Y plane (lbs) 

Aft lug reaction force X-direction (lbs) 

Forward lug reaction force X-direction (lbs) 

Aft lug reaction force Y-direction (lbs) 

Forward lug reaction force Y-direction (lbs) 

Aft lug reaction force Z-direction (lbs) 

Forward lug reaction force Z-direction (lbs) 
External moments X-direction PHIDD*IX (in-lbs) 
Combined aero and inertial moment roll-axis store 
coordinate system (in-lbs) 

External moments Y-direction THEDD*1IY (in=-lbe, 
Combined aero and inertial moment pitch-axis store 
coordinate system (in-lbs) 

External moments Z-direction PSIDD*1IZ (in-lbs) 
Combined aero and inertial moment yaw-axis store 
coordinate system (in-lbs) 

Acceleration X-direction (g) 

Acceleration Y-direction (9) 


Acceleration Z-direction (9) 
@ Angular acceleration X-axis (rad/sec**2) 


W Angular acceleration Z-axus (rad/sec**2) 
External loads X-direction W*NX (lbs) 


External airloads X-direction (lbs) 
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EXC 


EY 


PYA 


Eeuc 


PZA 


Sul 
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$3 


S4 


SS 


S6 


Combined airloads and inertial loads X-direction 
store coordinate system PX+PXA (lbs) 

External loads Y-direction W*NY 

External airloads Y-direction (lbs) 

Combined airloads and inertial loads Y-direction 
store coordinate system PY+PYA (lbs) 

External loads Z-direction W*NZ (lbs) 

External airloads Z-direction 

Combined airloads and inertial loads Z-direction 
Store coordinate system PZ+PZA (lbs) 

Dynamic pressure (lb/ft-sq) 

Reference radius (in) 

Distance between forward and aft sway brace pads 
(En) 

Distance store CG to fwd sway brace pad X-direction 
(na) 

Distance store CG to aft sway brace pad X-direction 
(in) 

Distance store CG to near sway brace pad 
YoGimectien (in) 

Distance store CG to far sway brace pad Y-direction 
(in). Note S3 equal to S4 when store symmetric 
Distance store CG to right sway brace pad 
Z-direction (in) 

Distance store CG to left sway brace pad 


Z-direction. Note S5 usually equal to S6 
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SUS 


3 LOLPAET 


SLZOTr WD 


SYART 


SYFWD 


SZAFT 


SZEFWD 


THEDD 


VEL 


W 


Reference area (ft-sq) 

Total aft sway brace load (lbs) 

Total forward sway brace load (lbs) 

Aft sway brace reaction force Y-direction (lbs) 
Forward sway brace reaction force Y-direction (lbs) 
Aft sway brace reaction force Z-direction (lbs) 


Forward sway brace reaction force Z-direction (lbs) 
6 Angular acceleration Y-axis (rad/sec**2) 


Aircraft velocity (ft/sec) 


Weight of store (lbs) 
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APPENDIX B: PROGRAM MIL8591H COMPUTER CODE* 


PROGRAM MIL8591 
C* PROGRAM TO COMPUTE STORE LOADS USING MIL-A-8591H 


C* DEFINITON OF TERMS 

C* W-WEIGHT OF STORE LBS 

C* IX-MASS MOMENT OF INERTIA X-AXIS SLUG-FT**2 

C* IY-MASS MOMENT OF INERTIA Y-AXIS SLUG-FT**2 

C* IZ-MASS MOMENT OF INERTIA Z-AXIS SLUG-FT**2 

C* ALL DISTANCES IN INCHES 

C* S- DISTANCE BETWEEN FORWARD AND AFT SWAY BRACE PADS 

C* S1-DISTANCE STORE CG TO FWD SWAY BRACE PAD X-DIRECTION 
C* S2-DISTANCE STORE CG TO AFT SWAY BRACE PAD X-DIRECTION 
C* S3-DISTANCE STORE CG TO NEAR SWAY BRACE PAD Y—-DIRECTION 
Cx S4-DISTANCE STORE CG TO FAR SWAY BRACE PAD Y-DIRECTION 
C* NOTE S3 EQUAL TO S4 WHEN STORE SYMMETRIC 

C* SS-DISTANCE STORE CG TO RIGHT SWAY BRACE PAD Z—-DIRECTION 
C* S6-DISTANCE STORE CG TO LEFT SWAY BRACE PAD Z—-DIRECTION 
C* NOTE S5 USUALLY EQUAL TO S6 

C* L=- DISTANCE BETWEEN LUGS 

C* L1-DISTANCE STORE CG TO FWD LUG X=-DIRECTION 

C* L2-DISTANCE STORE CG TO AFT LUG X-DIRECTION 

Cx L3-DISTANCE STORE CG TO FWD LUG Y-DIRECTION 

Cx L4-DISTANCE STORE CG TO AFT LUG Y-DIRECTION 

Cx NOTE L3 USUALLY EQUAL TO L4 

C* LS-DISTANCE STORE CG TO FWD LUG Z—-DIRECTION 

C* L6—-DISTANCE STORE CG TO AFT LUG Z-DIRECTION 

C* NOTE L5 USUALLY EQUAL TO L6 

Cx LXFWD-FORWARD LUG REACTION FORCE X-DIRECTION LBS 

C* LYFWD-FORWARD LUG REACTION FORCE Y-DIRECTION LBS 

C* LZFWD-FORWARD LUG REACTION FORCE Z-DIRECTION LBS 

Cx LXAFT-AFT LUG REACTION FORCE X-DIRECTION LBS 

C* LYAFT-AFT LUG REACTION FORCE Y-DIRECTION LBS 

C* LZAFT-AFT LUG REACTION FORCE Z-DIRECTION LBS 

Cx LTOTFWD- FORWARD LUG SHEAR FORCE X-Y PLANE LBS 

C* LTOTAFT—- AFT LUG SHEAR FORCE X-Y PLANE LBS 

Cx SYFWD- FORWARD SWAY BRACE REACTION FORCE Y-DIRECTION LBS 
C* SZFWD-FORWARD SWAY BRACE REACTION FORCE Z-DIRECTION LBS 
C* STOTFWD- TOTAL FORWARD SWAY BRACE LOAD LBS 

C* SYAFT-AFT SWAY BRACE REACTION FORCE Y-DIRECTION LBS 

C* SZAFT-AFT SWAY BRACE REACTION FORCE Z-DIRECTION LBS 

C* STOTAFT-— TOTAL AFT SWAY BRACE LOAD LBS 

C* NX-ACCELLERATION X—-DIRECTION 

C* NY-ACCELLERATION Y-DIRECTION 

C* NZ-ACCELLERATION Z—-DIRECTION 

C* PX-EXTERNAL LOADS X-DIRECTION W*NX 

C* PY-EXTERNAL LOADS Y-DIRECTION W*NY 


* 


copy of the computer program may be obtained by writing to 
the chairman of the Department of Aeronautics and Astronautics 
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C* PZ-EXTERNAL LOADS Z-DIRECTION W*NZ 
C* PHIDD-ANGULAR ACCELERATION X—-DIRECTION RAD/SEC**2 

C* THEDD-ANGULAR ACCELERATION Y-DIRECTION RAD/SEC**2 

C* PSIDD-ANGULAR ACCELERATION Z—-DIRECTION RAD/SEC**2 

C* MX-EXTERNAL MOMENTS X-DIRECTION PHIDD*I 

C* MY-EXTERNAL MOMENTS Y-DIRECTION THEDD*I 

C* MZ—-EXTERNAL MOMENTS Z—-DIRECTION PSIDD*I 

C* PXA-EXTERNAL AIRLOADS X-DIRECTION 

C* PYA-EXTERNAL AIRLOADS Y-DIRECTION 

C* PZA-EXTERNAL AIRLOADS Z—-DIRECTION 

C*  AMXA-AERODYNAMIC ROLL MOMENT 

C* AMYA-AERODYNAMIC PITCH MOMENT 

C* AMZA-AEBERODYNAMIC YAW MOMENT 

C* PXC-COMBINED AIRLOADS AND INERTIAL LOADS X-DIRECTION STORE COORD SYS 
C* PYC=-COMBINED AIRLOADS AND INERTIAL LOADS Y-DIRECTION STORE COORD SYS 
C* PZC-COMBINED AIRLOADS AND INERTIAL LOADS Z-DIRECTION STORE COORD SYS 
C* MXC-COMBINED AERO AND INERTIAL MOMENT ROLL-DIRECTION STORE COORD SYS 


C* MYC-COMBINED AERO AND INERTIAL MOMENT PITCH-DIRECTION STORE COORD SYS 
C* MZC-COMBINED AERO AND INERTIAL MOMENT YAW-DIRECTION STORE COORD SYS 


C* ANGA -— ANGLE OF ATTACK OF STORE (DEG) 

C* ANGS - ANGLE OF SIDE-SLIP (DEG) 

C* QUE - DYNAMIC PRESSURE (LB/FT-SQ) 

C* VEL - AIRCRAFT VELOCITY (FT/SEC) 

C* DEN - AIR DENSITY (SLUG/FT*%3) 

C* SR —- REFERENCE AREA (FT-SQ) 

C* R — REFERENCE RADIUS (IN) 

C* LR —- REFERENCE LENGTH (IN) 

C* Cl —- LIFT COEFFICIENT SLOPE (PER DEG) 

C* C2 — PITCH MOMENT COEFFICIENT SLOPE (PER DEG) 

C* C6 — SIDE FORCE COEFFICIENT SLOPE (PER DEG) 

C* C7 — YAW MOMENT COEFFICIENT SLOPE (PER DEG) 

C* CD —- DRAG COEFFICIENT 

C* EPS- ANGLE BETWEEN RADIUS OF CURVATURE OF STORE AND Z-AXIS (EPS) DEG. 
C* DEL1-ROTATION AROUND THE X-AXIS FOR STORES MOUNTED AT 
c* A ROLL ANGLE W.R.T. THE AIRCRAFT AXIS SYSTEM DEG 


C* THESE ARE VARIABLES USED IN THE SUBROUTINES TO CALCULATE LOADS 
INCLUDE ’ TOPJUNK.INC’ 
WRITE (6,*)’ENTER 1 FOR INTERACTIVE 2 FOR BATCH I0’ 
READ (5,*) NN 
OPEN (UNIT=11,FILE=’ FORO11.DAT’ , STATUS=’ UNKNOWN’ ) 
IF (NN.EQ.1) THEN 
NUMIN=5 
NUMOUT=6 
CALL GETDATA1 
ELSE 
NUMIN=10 
NUMOUT=11 
CALL GETDATA2 
ENDIF 
CALL AIRLOAD 
CALL ADDLOAD 
CALL ROTATE 
IF (NN.EQ.1) THEN 
NUMOUT=11 
ENDIF 
C* IF THE CG IS LATERALLY BETWEEN THE SWAY BRACES THEN PZ 
C* POSITIVE OR UP IS THE WORST GASE 
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C* IF THE CG OF THE STORE IS LATERALLY OUTSIDE THE SWAY BRACES 
C* A NEGATIVE OR DOWN PZ MUST ALSO BE ANALYSED 


br (GHes. EO 1 OR CHSE BOs FORe CHoe.bO.6) THEN 


IF (NX.LT.0.AND.NY.LT.0.AND.NZ.GT.0.AND.THEDD 
7 Gl .U.. AND. PSIDD-GT .O) THEN 
WRITE (NUMOUT, *) ’WORST CASE LOAD IS AT THE FWD LEFT SWAY BRACE’ 
PRINT *,’WORST CASE LOAD IS AT THE FWD LEFT SWAY BRACE’ 
CALL FSWYBRC 


ELSEIF (NX.LT.0.AND.NY.GT.0.AND.NZ.GT.0.AND. 

> THEDD.GT.0.AND.PSIDD.LT.O) THEN 
WRITE (NUMOUT, *) ’WORST CASE LOAD IS AT THE FWD RIGHT SWAY BRACE’ 
PRINT *,’WORST CASE LOAD IS AT THE FWD RIGHT SWAY BRACE’ 


CALL FSWYBRC 


ELSEIF (NX.GT.0.AND.NY.LT.0.AND.NZ.GT.0.AND.THEDD 
> .LT.0.AND.PSIDD.LT.O) THEN 


WRITE (NUMOUT, *) "WORST CASE LOAD IS AT THE AFT LEFT SWAY BRACE’ 
PRINT *,’WORST CASE LOAD IS AT THE AFT LEFT SWAY BRACE’ 
CALL ASWYBRC 


ELSEIF (NX.GT.0.AND.NY.GT.0.AND.NZ.GT.0.AND.THEDD 
> .LT.0.AND.PSIDD.GT.0O) THEN 


WRITE (NUMOUT, *) "WORST CASE LOAD IS AT THE AFT RIGHT SWAY BRACE’ 


PRINT *,’WORST CASE LOAD IS AT THE AFT RIGHT SWAY BRACE’ 
CALL ASWYBRC 


ELSEIF (NX.GT.0.AND.NZ.LT.0.AND.THEDD.LT.0) THEN 
WRITE (NUMOUT, *) ’WORST CASE LOAD IS AT THE FORWARD HOOK’ 
PRINT *,’WORST CASE LOAD IS AT THE FORWARD HOOK’ 

CALL FHOOK 


ELSEIF (NX.LT.0.AND.NZ.LT.0) THEN 

WRITE (NUMOUT, *) "WORST CASE LOAD IS AT THE AFT HOOK’ 
PRINT *,’WORST CASE LOAD IS AT THE AFT HOOK’ 

CALL AHOOK 


ELSE 
WRITE (NUMOUT, *) ’THE INERTIAL "G" FACTORS AND ANGULAR’ 
WRITE (NUMOUT, *) ‘ACCELERATIONS YOU ENTERED DO NOT PROVIDE’ 
WRITE (NUMOUT, *) ’A WORST CASE LOAD. CONSULT TABLE D-1;’ 
WRITE (NUMOUT, *) ‘APPENDIX D MIL-8591H.’ 
PRINT *,’THE INERTIAL "G" FACTORS AND ANGULAR’ 
PRINT *,’ACCELERATIONS YOU ENTERED DO NOT PROVIDE’ 
PRINT *,’A WORST CASE LOAD. CONSULT TABLE D-1;’ 
PRINT *,’ APPENDIX D MIL-8591H.’ 
ENDIF 


C* IF THE CG OF THE STORE IS LATERALLY OUTSIDE THE SWAY BRACES 
C* A NEGATIVE OR DOWN PZ MUST BE ANALYSED 


ELSEIF (CHSB.EQ.2.OR. CHSB.EQ.4.0OR. CHSB.EQ.5.OR.CHSB.EQ.7) THEN 
IF (NX.LT.0.AND.NY.LT.0.AND.NZ.GT.0.AND.PHIDD.LT.0.AND.THEDD 


Paeecn.0 AND. PSIDD. Gro) THEN 
WRITE (NUMOUT, *) ’WORST CASE LOAD IS AT THE FWD LEFT SWAY BRACE’ 
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PRINT *,’WORST CASE LOAD IS AT THE FWD LEFT SWAY BRACE’ 
CALL FSWYBRC 


ELSEIF (NX.LT.0.AND.NY .GT.0.AND.PHIPD 3G. 0 AND’: 

> THEDD.GT.0.AND.PSIDD.LT.O) THEN 
WRITE (NUMOUT, *) 72022 IO IOI I III II 
WRITE (NUMOUT, *) STORE CG IS LATERALLY OUTSIDE SWAY BRACES’ 
WRITE (NUMOUT, *) ‘THEREFORE AFT HOOK AND FWD SWAY BRACE’ 
WRITE (NUMOUT, *) “SHOULD BE ANALYSED’ 


WRITE (NUMOUT, *) KKK KK KK KKK KKK KKK KKK KKK KKK KKK KKK KK KKKKK 
PRINT *, % KRRKKKK KK KKK KKK KKK KKK KKK KKK KK KKK KK KKK KKKKKKKKKKKKKKS 


a 
PRINT *,’STORE CG IS LATERALLY OUTSIDE SWAY BRACES’ 
PRINT *,’THEREFORE AFT HOOK AND FWD SWAY BRACE’ 
PRINT *,’SHOULD BE ANALYSED’ 
PRINT cot KKK KKK KKK KKK KKK KKK KKK KK KK KKK KKK KKK KK KKK KKK KK KK KKK 
PRINT *,’ ¢ 
WRITE (NUMOUT, *) 
WRITE (NUMOUT, *) "WORST CASE LOAD IS AT THE AFT HOOK’ 
PRINT *,’WORST CASE LOAD IS AT THE AFT HOOK’ 
WRITE (NUMOUT, *) "WORST CASE LOAD IS AT THE FWD RIGHT SWAY BRACE’ 
PRINT *,’WORST CASE LOAD IS AT THE FWD RIGHT SWAY BRACE’ 
WRITE (NUMOUT, *) 
CALL FSWYBRC 
CALL AHOOK 


ELSEIF (NX.GT.0.AND.NY.LT.0.AND.PHIDD.LT.0.AND.THEDD 

> .LT.Q.AND.PSIDD.LT.O) THEN 
WRITE (NUMOUT, *) 122000802 233 IOI IOI III I III 
WRITE (NUMOUT, *)’STORE CG IS LATERALLY OUTSIDE SWAY BRACES’ 
WRITE (NUMOUT, *) ’ THEREFORE FWD HOOK AND AFT SWAY BRACE’ 
WRITE (NUMOUT, *) ‘SHOULD BE ANALYSED’ 


WRITE (NUMOUT, *) GKKKKK KKK KKK KKK KKK KK KKK KKK KK KKK KKK KKKKKKKK KKK 
PRINT * KKK KKK KKK KKK KK KKK KKK KKK KKK KKK KK KKK KKK KKKKKKKKKKKK 


tA 
PRINT *,’STORE CG IS LATERALLY OUTSIDE SWAY BRACES’ 
PRINT *,’THEREFORE FWD HOOK AND AFT SWAY BRACE’ 
PRINT *,’SHOULD BE ANALYSED’ 
PRINT oi G—KKKKKKKKKKKKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKKS 
WRITE (NUMOUT, *) 
PRINT *,/ ¢ 
WRITE (NUMOUT, *) "WORST CASE LOAD IS AT THE AFT LEFT SWAY BRACE’ 
PRINT *,’WORST CASE LOAD IS AT THE AFT LEFT SWAY BRACE’ 
WRITE (NUMOUT, *) ‘WORST CASE LOAD IS AT THE FORWARD HOOK’ 
PRINT *,’WORST CASE LOAD IS AT THE FORWARD HOOK’ 
CALL ASWYBRC 
CALL FHOOK 


ELSEIF (NX.GT.0.AND.NY.GT.0.AND.NZ.GT.0.AND.PHIDD.GT.0.AND.THEDD 
> .LT.0.AND.PSIDD.GT.0O) THEN 

WRITE (NUMOUT, *) WORST CASE LOAD IS AT THE AFT RIGHT SWAY BRACE’ 

PRINT *,’WORST CASE LOAD IS AT THE AFT RIGHT SWAY BRACE’ 

CALL ASWYBRC 


ELSEIF (NX.GT.0.AND.NZ.LT.0.AND.THEDD.LT.0) THEN 

WRITE (NUMOUT, *) “WORST CASE LOAD IS AT THE FORWARD HOOK’ 
PRINT *,’WORST CASE LOAD IS AT THE FORWARD HOOK’ 

CALL FHOOK 


ELSEIF (NX.LT.0.AND.NZ.LT.0) THEN 

WRITE (NUMOUT, *) ‘WORST CASE LOAD IS AT THE AFT HOOK’ 
PRINT *,’WORST CASE LOAD IS AT THE AFT HOOK’ 

CALL AHOOK 
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ELSE 

WRITE (NUMOUT, *) ’THE INERTIAL "G" FACTORS AND ANGULAR’ 
WRITE (NUMOUT, *) ’ ACCELERATIONS YOU ENTERED DO NOT PROVIDE’ 
WRITE (NUMOUT, *)’A WORST CASE LOAD. CONSULT TABLE D-1;’ 
WRITE (NUMOUT, *) “APPENDIX D MIL-8591H.’ 

PRINT *,’THE INERTIAL "G" FACTORS AND ANGULAR’ 


a 
PRINT *,’ ACCELERATIONS YOU ENTERED DO NOT PROVIDE’ 
PRINT *,’A WORST CASE LOAD. CONSULT TABLE D-1;’ 
PRINT *,’ APPENDIX D MIL—-8591H.’ 
ENDIF 
ENDIF 


CALL OUTPUT 


oe 


IF (NN.EQ.1) THEN 
NUMOUT=6 

CALL OUTPUT 
ENDIF 

IF (NN.EQ.2) THEN 
NUMOUT=6 

CALL OUTPUT 
ENDIF 

SLOP 

END 


SUBROUTINE OUTPUT 
INCLUDE ’ TOPJUNK.INC’ 
WRITE (NUMOUT, *) 


WRITE (NUMOUT, *) ’ **TNERT PAI LOADS TAPPLIED*=* 
WRITE (NUMOUT, *) ” LOAD FACTOR G"S ANGULAR ACCELERATION 
RAD/SEC**2’ 

WRITE (NUMOUT, *) ’ LONGITUDINAL=’ ,NX,’ ROLL=’ , PHIDD 

WRITE (NUMOUT, *) ’ LATERAL mY, © fleH—" THEDD 

WRITE (NUMOUT, *) ‘ VERTICAL =',NZ,' YAW=’,PSIDD 


WRITE (NUMOUT, *) 
WRITE (NUMOUT, *) 


WRITE (NUMOUT, *) ’ **AERODYNAMIC LOADS APPLIED**’ 
WRITE (NUMOUT, *) ’ AIRLOAD LB AERODYNAMIC MOMENT 
IN-LB’ 


’,PXA,’ ROLL =’, AMXA 
’, PYA,’ PITCH =’ , AMYA 
f, PZA,’ YAW =',AMZA 


WRITE (NUMOUT, *) ’ LONGITUDINAL 
WRITE (NUMOUT, *) ’ LATERAL 
WRITE (NUMOUT, *) ‘ VERTICAL 
WRITE (NUMOUT, *) 

WRITE (NUMOUT, *) 

WRITE (NUMOUT, *) STORE WEIGHT LBS =’',W 
WRITE (NUMOUT, *) “MASS MOMENT OF INERTIA X-AXIS LB-IN-SEC*2 =',1X 
WRITE (NUMOUT, *) "MASS MOMENT OF INERTIA Y-AXIS LB-IN-SEC*2 =’, IY 
WRITE (NUMOUT, *) "MASS MOMENT OF INERTIA Z-AXIS LB-IN-SEC*2 =',12Z 
WRITE (NUMOUT, *) 

WRITE (NUMOUT, *) 

WRITE (NUMOUT, *) "HOOK/LUG SPACING IN. =’,L 

WRITE (NUMOUT, *) THE STORE CG IS’,Li,’IN FROM THE FWD HOOK/LUG’ 
WRITE (NUMOUT, *) 

WRITE (NUMOUT, *) “LONGITUDINAL SWAY BRACE SPACING IN =’,S 

WRITE (NUMOUT, *) THE STORE CG IS’,S1,’IN FROM THE FWD SWAY BRACE’ 
(CHAL.EQ. 3) THEN 

WRITE (NUMOUT, *) 

WRITE (NUMOUT, *) "AIRLOADS WERE CALCULATED CASE 3’ 
WRITE (NUMOUT, *) 

WRITE (NUMOUT, *) “AIRCRAFT VELOCITY FT/SEC =", VEL 
WRITE (NUMOUT, *) "STORE ANGLE OF ATTACK DEG =’,ANGA 
WRITE (NUMOUT, *) "STORE ANGLE OF SIDESLIP DEG =’,ANGS 
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ENDIE 


RETURN 
END 


SUBROUTINE GETDATA1 
INCLUDE ’ TOPJUNK.INC’ 


WRITE (NUMOUT,*)’ INPUT STORE CG POSITION WITH RESPECT TO 
> SWAY BRACE PADS. ENTER THE NUMBER 


> CORRESPONDING TO THE CONDITION’ 
WRITE (NUMOUT. =). 4 — STORE CG BETWEEN SWAY BRACE PADS 
LONGITUDINALLY 


> AND LATERALLY’ 


WRITE (NUMOUT, *) 


WRITE (NUMOUT,*) 2— STORE CG BETWEEN SWAY BRACE PADS 
LONGITUDINALLY 
> BUT OUTSIDE SWAY BRACE PADS LATERALLY’ 


Vv 


WRITE (NUMOUT, *) 
WRITE (NUMOUT, *)’3- STORE CG FORWARD OF FRONT SWAY BRACE PAD 


LONGITUDINALLY, BUT BETWEEN SWAY BRACE PADS LATERALLY’ 


WRITE (NUMOUT, *) 
WRITE (NUMOUT, *)’4—- STORE CG FORWARD OF FRONT SWAY BRACE PAD 


> LONGITUDINALLY, AND OUITSIDE SWAY BRACE PADS LATERALLY’ 
WRITE (NUMOUT, *) 
WRITE (NUMOUT,*)’5—- STORE CG BETWEEN SWAY BRACE PAD 

> LONGITUDINALLY, BUT OUTSIDE SWAY BRACE PADS LATERALLY’ 
WRITE (NUMOUT, *) 
WRITE (NUMOUT, *)’6—- STORE CG AFT OF REAR SWAY BRACE PAD 

> LONGITUDINALLY, BUT BETWEEN SWAY BRACE PADS LATERALLY’ 
WRITE (NUMOUT, *) 
WRITE (NUMOUT, *)’7— STORE CG AFT OF REAR SWAY BRACE PAD 

> LONGITUDINALLY, AND OUTSIDE SWAY BRACE PADS LATERALLY’ 


READ (NUMIN, *) CHSB 
IF (CHSB_EQaj)) THEN 


WRITE (NUMOUT, *)’i-— STORE CG BETWEEN SWAY BRACE PADS 


LONGITUDINALLY 


> 


AND LATERALLY’ 

WRITE (NUMOUT, *) 

ELSE IE GHoB. EO22) THEN 

WRITE (NUMOUT, *) ’ 2- STORE CG BETWEEN SWAY BRACE PADS 


LONGITUDINALLY 


> 


BUT OUTSIDE SWAY BRACE PADS LATERALLY’ 


WRITE (NUMOUT, *) 
ELSEIF (CHSB.EQ.3) THEN 
WRITE (NUMOUT, *)’3- STORE CG FORWARD OF FRONT SWAY BRACE PAD 


LONGITUDINALLY, BUT BETWEEN SWAY BRACE PADS LATERALLY’ 

WRITE (NUMOUT, *) 

ELSEIF (CHSB.EQ.4) THEN 

WRITE (NUMOUT, *)’4— STORE CG FORWARD OF FRONT SWAY BRACE PAD 
LONGITUDINALLY, AND OUITSIDE SWAY BRACE PADS LATERALLY’ 
WRITE (NUMOUT, *) 

ELSEIF (CHSB.EQ.5) THEN 


WRITE (NUMOUT, *)’5—- STORE CG BETWEEN SWAY BRACE PAD 
LONGITUDINALLY, BUT OUTSIDE SWAY BRACE PADS LATERALLY’ 
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Cc* 


Cc* 


Cx 


WRITE (NUMOUT, *) 
ELSEIF (CHSB.EQ.6) THEN 
WRITE (NUMOUT,*)’6— STORE CG AFT OF REAR SWAY BRACE PAD 


> LONGITUDINALLY, AND BETWEEN SWAY BRACE PADS LATERALLY’ 


WRITE (NUMOUT, *) 

ELSEIF (CHSB.EQ.7) THEN 

WRITE (NUMOUT,*)’'7— STORE CG AFT OF REAR SWAY BRACE PAD 
> LONGITUDINALLY, AND OUTSIDE SWAY BRACE PADS LATERALLY’ 


ENDIF 


WRITE (NUMOUT, *)’ INPUT STORE CG POSITION WITH RESPECT TO 
>HOOK/LUGS. ENTER THE NUMBER CORRESPONDING TO THE CONDITION’ 

WRITE (NUMOUT, *)’1— STORE CG BETWEEN HOOK/LUGS’ 

WRITE (NUMOUT, *)’2-—- STORE CG IN FRONT OF FORWARD HOOK/LUGS’ 

WRITE (NUMOUT, *)’3— STORE CG AFT OF THE REAR HOOK/LUGS’ 


READ (NUMIN,*) CHHL 
WRITE (NUMOUT, *) 
IF (CHHL.EQ.1) THEN 
WRITE (NUMOUT, *)’1—- STORE CG BETWEEN HOOK/LUGS’ 
ELSEIF (CHHL.EQ.2) THEN 
WRITE (NUMOUT, *)’2— STORE CG IN FRONT OF FORWARD HOOK/LUGS’ 
ELSEIF (CHHL.EQ.3) THEN 
WRITE (NUMOUT, *)’3- STORE CG AFT OF REAR HOOK/LUGS’ 
ENDIF 
WRITE (NUMOUT, *) 
WRITE (NUMOUT, *) 
WRITE (NUMOUT,*)’ INPUT STORE WEIGHT (W) LBS’ 
READ (NUMIN,*) W 
WRITE (NUMOUT, *)’ INPUT STORE MOMENT OF INERTIA 
> X-AXIS (IX) SLUG-FT**2! 
READ (NUMIN,*) IX 
CONVERT TO LB-IN-SEC**2 
IX=12* 1X 
WRITE (NUMOUT,*) INPUT STORE MOMENT OF INERTIA 
> Y-AXIS (IY) SLUG-FT**2’ 
READ (NUMIN,*) IY 


CONVERT TO LB-IN-SEC**2 

7 ay, 
WRITE (NUMOUT, *)’ INPUT STORE MOMENT OF INERTIA 
pacer bo (127)  ShLUG-PT**2" 
READ (NUMIN,*) IZ 


CONVERT TO LB-IN-SEC**2 

IZ=12*1Z 
WRITE (NUMOUT,*) ‘INPUT DISTANCE BETWEEN FOR AND 
> AFT SWAY BRACE PADS (S) IN’ 
READ (NUMIN,*) S 
WRITE (NUMOUT, *)’ INPUT DISTANCE BETWEEN CG AND 
> FOR SWAY BRACE PADS X=DIRECTION (S1) IN’ 
READ (NUMIN,*) Sl 
WRITE (NUMOUT,*)’ INPUT DISTANCE BETWEEN CG AND 
> AFT SWAY BRACE PADS X-DIRECTION (S2) IN’ 
READ (NUMIN,*) S2 
WRITE (NUMOUT,*)’ INPUT DISTANCE BETWEEN CG AND 
> NEAR SWAY BRACE PADS Y-DIRECTION (S3) IN’ 
READ (NUMIN,*) S3 
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WRITE (NUMOUT,*)’ INPUT DISTANCE BETWEEN CG AND 
> FAR SWAY BRACE PADS Y-DIRECTION (S4) IN’ 
READ (NUMIN,*) $4 
WRITE (NUMOUT,*)’ INPUT DISTANCE BETWEEN CG AND 
> RIGHT SWAY BRACE PADS Z-DIRECTION (S5) IN’ 
READ (NUMIN,*) $5 
WRITE (NUMOUT,*)’ INPUT DISTANCE BETWEEN CG AND 
> LEFT SWAY BRACE PADS Z-DIRECTION (S6) IN’ 
READ (NUMIN,*) S6 
WRITE (NUMOUT, *)’ INPUT DISTANCE BETWEEN FOR AND 
> AFT LUGS X-DIRECTION (L) IN’ 
READ (NUMIN,*) L 
WRITE (NUMOUT, *)’ INPUT DISTANCE BETWEEN CG AND 
> FWD LUG X-DIRECTION (L1) IN’ 
READ (NUMIN,*) Ll 
WRITE (NUMOUT, *)’ INPUT DISTANCE BETWEEN CG AND 
> AFT LUG X-DIRECTION (L2) IN’ 
READ (NUMIN,*) L2 
WRITE (NUMOUT,*)’ INPUT DISTANCE BETWEEN CG AND 
> FWD LUG Y-DIRECTION (L3) IN’ 
READ (NUMIN,*) L3 
WRITE (NUMOUT, *) “INPUT DISTANCE BETWEEN CG AND 
> AFT LUG Y-DIRECTION (L4) IN’ 
READ (NUMIN,*) L4 
WRITE (NUMOUT, *) INPUT DISTANCE BETWEEN CG AND 
> FWD LUG Z-DIRECTION (L5) IN’ 
READ (NUMIN,*) L5 
WRITE (NUMOUT,*)’ INPUT DISTANCE BETWEEN CG AND 
> AFT LUG Z-DIRECTION (L6) IN’ 
READ (NUMIN,*) L6 
WRITE (NUMOUT, *)’ INPUT EXTERNAL LOAD FACTOR X-DIRECTION NX’ 
READ (NUMIN,*) NX 
WRITE (NUMOUT,*)’ INPUT EXTERNAL LOAD FACTOR Y-DIRECTION NY’ 
READ (NUMIN,*) NY 
WRITE (NUMOUT, *)’ INPUT EXTERNAL LOAD FACTOR Z-DIRECTION NZ’ 
READ (NUMIN,*) NZ 
WRITE (NUMOUT,*)’ INPUT EXTERNAL ANGULAR ACCELLERATION 
> X=-DIRECTION (PHIDD) RAD/SEC*2’ 
READ (NUMIN,*) PHIDD 
WRITE (NUMOUT, *)’ INPUT EXTERNAL ANGULAR ACCELLERATION 
> Y-DIRECTION (THEDD) RAD/SEC*2’ 
READ (NUMIN,*) THEDD 
WRITE (NUMOUT, *)’ INPUT EXTERNAL ANGULAR ACCELLERATION 
> Z-DIRECTION (PSIDD) RAD/SEC%*2’ 
READ (NUMIN,*) PSIDD 
WRITE (NUMOUT, *) ‘INPUT ANGLE BETWEEN RADIUS OF CURVATURE 
> OF STORE AND Z-AXIS (EPS) DEG’ 
READ (NUMIN,*) EPS 
WRITE (NUMOUT,*) ‘INPUT THE ROLL ANGLE BETWEEN THE AIRCRAFT AND 
> THE STORE (DEL1) DEGREES’ 
READ (NUMIN, *) DEL1 


CONVERT DEGREES TO RADIANS 
EPS=EPS/57.2958 
DEL1=DEL1/57.2958 


H=(l5=8 9) ~TANICBES) 
K=H+ (L3+S3) 


WRITE (NUMOUT, *)’ INPUT HOW AIRLOADS ARE TO BE ANALYSED. 
> ENTER THE NUMBER CORRESPONDING TO THE METHOD DESIRED.’ 
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WRITE (NUMOUT, *)’1—-NEGLECT AIRLOADS (IE. CATS TRAPS)’ 
WRITE (NUMOUT,*)’2-INPUT AIRLOADS DIRECTLY FROM 

> WIND TUNNEL DATA’ 
WRITE (NUMOUT, *)’3-CALCULATE AIRLOADS’ 


READ (NUMIN,*) CHAL 


IF (CHAL.EQ.1) THEN 
WRITE (NUMOUT, *)’1—-NEGLECT AIRLOADS (IE. CATS TRAPS) ’ 
ELSEIF (CHAL.EQ.2) THEN 
WRITE (NUMOUT, *)’2-INPUT AIRLOADS DIRECTLY FROM 
> WIND TUNNEL DATA’ 
ELSEIF (CHAL.EQ.3) THEN 
WRITE (NUMOUT, *) ’3-CALCULATE AIRLOADS’ 
ENDIF 
RETURN 
END 


SUBROUTINE GETDATA2 

INCLUDE ’ TOPJUNK.INC’ 
READ (NUMIN, ’ (A)’) COM 
READ (NUMIN,*) CHSB 

READ (NUMIN, ’ (A)’) COM 

READ (NUMIN,*) CHHL 

READ (NUMIN,’ (A)’) COM 
READ (NUMIN,*) W 

READ (NUMIN, ’ (A)’) COM 
READ (NUMIN,*) IX 

READ (NUMIN,’ (A)’) COM 
READ (NUMIN,*) IY 


READ (NUMIN, ’ (A)’) COM 
READ (NUMIN,*) IZ 


READ (NUMIN,’ (A)’) COM 
READ (NUMIN,*) S 

READ (NUMIN, ’ (A)’) COM 
READ (NUMIN,*) $1 

READ (NUMIN, ’ (A)’) COM 
READ (NUMIN,*) S2 

READ (NUMIN, ’ (A)’) COM 
READ (NUMIN,*) $3 

READ (NUMIN,’ (A)’) COM 
READ (NUMIN,*) S4 

READ (NUMIN, ’ (A)’) COM 
READ (NUMIN,*) S5 

READ (NUMIN, ’ (A)’) COM 
READ (NUMIN,*) S6 

READ (NUMIN, ’ (A)’) COM 
READ (NUMIN,*) L 

READ (NUMIN,’ (A)’) COM 
READ (NUMIN,*) Ll 

READ (NUMIN,’ (A)’) COM 
READ (NUMIN,*) L2 

READ (NUMIN,’ (A)’) COM 
READ (NUMIN,*) L3 

READ (NUMIN,’ (A)’) COM 
READ (NUMIN,*) L4 

READ (NUMIN,’ (A)’) COM 

READ (NUMIN,*) L5 

READ (NUMIN,’ (A)’) COM 


a 


READ (NUMIN,*) L6 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN,*) NX 
READ (NUMIN, ’ (A)’) COM 
READ (NUMIN,*) NY 
READ (NUMIN, ’ (A)’) COM 
READ (NUMIN,*) NZ 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN,*) PHIDD 
READ (NUMIN, ’ (A)’) COM 
READ (NUMIN,*) THEDD 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN,*) PSIDD 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN,*) EPS 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN, *) DEL1 
READ (NUMIN, ’ (A)’) COM 
READ (NUMIN,*) CHAL 
READ (NUMIN, ’ (A) ’) COM 
READ (NUMIN, ’ (A)’) COM 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN, ’ (A)’) COM 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN, ’ (A)’) COM 
IF (CHAL.EQ.2) THEN 
READ (NUMIN, *) PXA 
READ (NUMIN, ’ (A)’) COM 
READ (NUMIN, *) PYA 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN, *)PZA 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN, *) AMXA 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN, *) AMYA 
READ (NUMIN, ’ (A)’) COM 
READ (NUMIN, *) AMZA 
ENDIF 


BYPASSES READING DIRECT INPUT AIRLOADS 
IF (CHAL.EQ.3) THEN 
DO I=1,6 

READ (NUMIN, *) PXA 

READ (NUMIN, ’ (A) ’) COM 
ENDDO 


READ (NUMIN, ’ (A)’) COM 
READ (NUMIN, ’ (A)’) COM 
READ (NUMIN, ‘ (A) ’) COM 
READ (NUMIN, *) VEL 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN, *) ANGA 
READ (NUMIN, ’ (A)’) COM 
READ (NUMIN, *) ANGS 
READ (NUMIN, ’ (A)’) COM 
READ (NUMIN, *) DEN 
READ (NUMIN, ’ (A) ’) COM 
READ (NUMIN, *)R 

READ (NUMIN, ’ (A)’) COM 
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READ (NUMIN,’ (A)’) COM 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN, *)SR 
READ (NUMIN, ’ (A)’) COM 
READ (NUMIN, *) LR 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN, *)Cl 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN, *)C2 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN, *)C6 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN, *)C7 
READ (NUMIN,’ (A)’) COM 
READ (NUMIN, *) CD 


ENDIF 
CONVERT IX FROM SLUG-FT**2 TO LB-IN-SEC**2 
ox—12* 07x 
TY=12*IyY 
7 —1 2 * 12, 

CONVERT DEGREES TO RADIANS 
EPS=EPS/57.2958 
DEL1=DEL1/57.2958 

H= (L5-S5) *TAN (EPS) 

K=H+ (L3+S3) 

RETURN 
END 


SUBROUTINE CACULATES AIRLOADS TO BE ANALYSED 
SUBROUTINE AIRLOAD 
INCLUDE 'TOPJUNK.INC’ 
FOR A BATCH FILE INPUTS ARE IN SUBR GETDATA2 
IF (NN.EQ.1) THEN 
IF (CHAL.EQ.1) THEN 
PXA=0 
PYA=0 
PZA=0 
AMXA=0 
AMYA=0 
AMZA=0 
ELSEIF (CHAL.EQ.2) THEN 
WRITE (NUMOUT, *)’ INPUT AIRLOADS DIRECTLY FROM WIND TUNNEL DATA’ 
WRITE (NUMOUT,*)’ INPUT LONGITUDINAL AIRLOAD (PXA) 
> LB X-DIRECTION POSITIVE AFT’ 
READ (NUMIN,*) PXA 
WRITE (NUMOUT,*)’ INPUT LATERAL AIRLOAD (PYA) 
> LB Y-DIRECTION POSITIVE LEFT’ 
READ (NUMIN,*) PYA 
WRITE (NUMOUT, *) INPUT VERTICAL AIRLOAD (PYZ) 
> LB Z-DIRECTION POSITIVE UP’ 
READ (NUMIN,*) PZA 
WRITE (NUMOUT, *)’ INPUT LATERAL MOMENT= ROLL (AMXA) 
> IN-LB POSITIVE CCW LOOKING FORWARD’ 
READ (NUMIN,*) AMXA 
WRITE (NUMOUT,*)’ INPUT LONGITUDINAL MOMENT— PITCH (AMYA) 
> IN-LB POSITIVE NOSE UP’ 
READ (NUMIN,*) AMYA 
WRITE (NUMOUT,*)’ INPUT VERTICAL MOMENT— YAW (AMZA) 
> IN-LB POSITIVE NOSE LEFT’ 
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READ (NUMIN, *) AMZA 
ELSEIF (CHAL.EQ.3) THEN 


WRITE (NUMOUT, *)’CALCULATE AIRLOADS FROM AERODYNAMIC DATA’ 


DEFINITION OF TERMS 
ANGA — ANGLE OF ATTACK OF STORE (DEG) 


ANGS - ANGLE OF SIDE-SLIP (DEG) 
QUE - DYNAMIC PRESSURE (LB/FT-SQ) 
VEL -— AIRCRAFT VELOCITY (FT/SEC) 
DEN — AIR DENSITY (SLUG/FT*3) 

SR - REFERENCE AREA (FT-SQ) 

LR — REFERENCE LENGTH (IN) 
R-REFERENCE RADIUS (IN) 


Cl = LIFT COnPPICIENT SLOPE (PER DEG) 


C2 PITCH MOMENT COEFFICIENT SLOPE (PER DEG) 


C6 —- SIDE FORCE COEFFICIENT SLOPE (PER DEG) 
C7 — YAW MOMENT COEFFICIENT SLOPE (PER DEG) 


58, DRAG COEFFICIENT 


WRITE (NUMOUT, *) 

WRITE (NUMOUT,*) ‘’ENTER THE AIRCRAFT CARRIAGE CONDITIONS 
> AND ATMOSPHERIC DATA AS PROMPTED’ 

WRITE (NUMOUT, *) 


WRITE (NUMOUT, *) ’AIRCRAFT VELOCITY - FT/SEC ’ 
READ (NUMIN, *) VEL 

WRITE (NUMOUT, *) 

WRITE (NUMOUT,*)’ ANGLE OF ATTACK -— DEG’ 

READ (NUMIN, *) ANGA 

WRITE (NUMOUT, *) 

WRITE (NUMOUT,*)’ANGLE OF SIDE-SLIP -— DEG’ 
READ (NUMIN, *) ANGS 

WRITE (NUMOUT, *) 

WRITE (NUMOUT,*)’AIR DENSITY -— SLUG/CUBIC FT ’ 
READ (NUMIN, *) DEN 

WRITE (NUMOUT, *)’REFERENCE RADIUS - IN’ 

READ (NUMIN, *)R 


WRITE (NUMOUT, *) 

WRITE (NUMOUT, *) ‘WRITE STORE AERODYNAMIC CHARACTERISTICS 
> (FOR UNKNOWN VALUES, ENTER A "0" AND DEFAULT VALUES 
> BASED ON A CYLINDRICAL SHAPE WILL BE ENTERED’ 

WRITE (NUMOUT, *) 

WRITE (NUMOUT, *) REFERENCE AREA —- FT-SQ’ 

READ (NUMIN, *) SR 

WRITE (NUMOUT, *) 


WRITE (NUMOUT, *)’REFERENCE LENGTH - IN’ 

READ (NUMIN, *) LR 

WRITE (NUMOUT, *) 

WRITE (NUMOUT, *)’ LIFT COEFFICIENT — PER DEG’ 
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Cx 


Cx 
Cx 


Cx 


Cx 


Cx 


READ (NUMIN, *)C1 
WRITE (NUMOUT, *) 


WRITE (NUMOUT,*)* PITCH MOMENT “COEFFICIENT = PER DEG’ 


READ (NUMIN, *) C2 
WRITE (NUMOUT, *) 


WRITE (NUMOUT,*)* SIDE—FORCE COEFFICIENT» = PER DEG" 


READ (NUMIN, *) C6 
WRITE (NUMOUT, *) 


WRITE (NUMOUT, *)’ YAW MOMENT COEFFICIENT - PER DEG’ 


READ (NUMIN, *) C7 
WRITE (NUMOUT, *) 


WRITE (NUMOUT, *) ‘DRAG COEFFICIENT ’ 


READ (NUMIN, *) CD 
ENDIF 
ENDIF 


DETERMINE DEFAULT VALUES 


IF (CHAL.EQ.3) THEN 


IF (SR.EQ.0) THEN 


SR = 3.14159 * ((2 * R / 12)** 2.) / 4 


ENDIF 
IF (LR.EQ.0) THEN 
LR = 22 * R 
ENDIF: 
“REM L = 11*D 
FROM HOERNER P.19-8 
IF (C1.EQ.0) THEN 
Gl = -.032 


‘REM HOERNER LIFT P. 


IF (C2.EQ.0) THEN 
25> .013 


‘REM HOERNER LIFT P. 


Tb -s(C6.EO-0) THEN 
€6.= Cl 


ENDIF 
Tene? 50> 0)o LHEN 
C7 =.C2 

ENDIF 


[TF (CD.EO.0) THEN 
CD = .3 


19-8 


1S ie: 


‘REM HOERNER DRAG P. 3-12 FIG. 21 


ENDIF 


QUE = DEN * VEL ** 2. / 2 


paaA = CD * OUE * SR 
‘REM DRAG FORCE 
PYA = Co * ANGS * QUE * SR 
‘REM SIDE FORCE ACTING AT .57*LR ASSUME CG (FROM HOERNER) 
RPZA = Cl * ANGA * QUE * SR 
‘REM LIFT FORCE ACTING AT .57*LR ASSUME CG (FROM HOERNER) 
AMXA = 0 
‘REM ROLL MOMENT 
AMYA = C2 * ANGA * QUE * SR * LR 
‘REM PITCH MOMENT 
AMZA = C/ * ANGS * QUE * SR 49LR 
‘REM YAW MOMENT 
ENDIF 


SUBROUTINE TO ADD INERTIAL LOADS AND AIRLOADS 
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SUBROUTINE ADDLOAD 
INCLUDE ’ TOPJUNK.INC’ 
PX=ABS (NX*W+P XA) 
PY=ABS (NY *W+PYA) 
PZ=ABS (NZ*W+P ZA) 
MX=ABS (PHIDD* IX+AMXA) 
MY=ABS (THEDD* IY+AMYA) 
MZ=ABS (PSIDD*1IZ+AMZA) 
RETURN 

END 


C* SUBROUTINE TO RESOLVE STORE COORDINATE SYSTEM WITH AIRCRAFT 
C* COORDINATE SYSTEM 
SUBROUTINE ROTATE 
INCLUDE ’ TOPJUNK.INC’ 
IF (DEL1.EQ.0) THEN 
EAC=E XA 
PYC=PY 
PZC-E2Z 
MXC=MX 
MYC=MY 
MZC=MZ 
ELSE 
PXC=PX 
PYC=PY*COS (DEL1) +PZ*SIN (DEL1) 
PZC=PZ*COS (DEL1) -PY*SIN (DEL1) 
MXC=MX 
MYC=MY*COS (DEL1) +MZ*SIN (DEL1) 
MZC=MZ*COS (DEL1) —-MY*SIN (DEL1) 
ENDIF 
RETURN 
END 


C* SUBROUTINE CALCULATES FORWARD SWAY BRACE LOADS 
SUBROUTINE FSWYBRC 
INCLUDE ’TOPJUNK.INC’ 
A=L5/ (2* (S1+L2) )+L3/ (S*TAN (EPS) ) 
B=S2*L5/ (S*K) 
C= (S275) +bs7/ K) 


IF (CHSB.EQ.2) THEN 

C= (L2*L3) / (L* (H+S4=L3) ) 
ELSEIF (CHSB.EQ.3) THEN 
C=(L2/ (L2-S1) ) * (.5+L3/K) 
ELSEIF (CHSB.EQ.4) THEN 

C= (L2*L3) / ((L2-S1) * (H+S4-L3) ) 
ENDIF 


D=S2/ (S*K) 
E=.5/(S1+1L2) 
F=1/(S* (TAN(EPS) ) ) 
SZFWD=A*PXC+B*PYC+C*PZC+D*MXC+E*MYC+FE *MZC 
SYFWD=SZFWD*TAN (EPS) 
STOTFWD=SZFWD/ (COS (EPS) ) 
WRITE (NUMOUT, *) 
WRITE (NUMOUT, *)’ SWAY BRACE LOAD LBS’ 
WRITE (NUMOUT, *) ‘’ VERTICAL=’ , SZFWD 
WRITE (NUMOUT, *) ’ LATERAL =’, SYFWD 
WRITE (NUMOUT, *) ’ TOTAL =', STOTFWD 
NUMOUT=6 
WRITE (NUMOUT, *) 
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WRITE (NUMOUT, *)’ SWAY BRACE LOAD LBS’ 
WRITE (NUMOUT, *) ‘’ VERTICAL=’ , SZFWD 
WRITE (NUMOUT, *) ‘LATERAL ='’,SYFWD 
WRITE (NUMOUT, *) ’ TOTAL =", STOTEWD 

NUMOUT=11 

RETURN 

END 


C* SUBROUTINE CALCULATES AFT SWAY BRACE LOADS 
SUBROUTINE ASWYBRC 
INCLUDE ’TOPJUNK.INC’ 
A=L5/ (2* (S2+L1) ) +L3/ (S*TAN (EPS) ) 
Bool Woy (s7K) 
C=(S1/S) * (.5+L3/K) 


IF (CHSB.EQ.5) THEN 
C= (L1*L3) / (L* (H+S4-—L3) ) 
ELSEIF (CHSB.EQ.6) THEN 
C=(L1* (L1—-S2) ) * (.5+L3/K) 
ELSEIF (CHSB.EQ.7) THEN 
C=(L1*L3) / ((L1—-S2) * (H+S4-L3) ) 
ENDIF 


D=S1/ (S*K) 

E=.5/(S2+L1) 

F=1/(S* (TAN (EPS) ) ) 
SZAFT=A*PXC+B*PYC+C*PZC+D *MXC+E*MYC+E *MZC 
SYAFT=SZAFT*TAN (EPS) 
STOTAFT=SZAFT/ (COS (EPS) ) 

WRITE (NUMOUT, *) 

WRITE (NUMOUT, *)’ SWAY BRACE LOAD LBS’ 
WRITE (NUMOUT, *) ‘ VERTICAL=’ , SZAFT 
WRITE (NUMOUT, *) LATERAL =’,SYAFT 
WRITE (NUMOUT, *) ‘TOTAL =', STOTAFT 
NUMOUT=6 

WRITE (NUMOUT, *) 

WRITE (NUMOUT, *)’ SWAY BRACE LOAD LBS’ 
WRITE (NUMOUT, *) ‘VERTICAL=’ , SZAFT 
WRITE (NUMOUT, *) “LATERAL =’,SYAFT 
WRITE (NUMOUT, *) ‘ TOTAL =', STOTAFT 
NUMOUT=11 

RETURN 

END 


C* SUBROTINE TO CALCULATE FORWARD HOOK LOAD 
SUBROUTINE FHOOK 
INCLUDE ’ TOPJUNK.INC’ 
AY=0 
BY=(S2/S) * ((L5* (TAN (EPS) ) /K)-1) 
C* FIND SIGN OF BY 
T=S5*TAN (EPS) — ($3+L3) 
IF (I.LT.0) THEN 
BY=-BY 
ELSEIF (I.GT.0) THEN 
BY=BY 
ELSE 
BY=0 
ENDIF 
CY=L2*L3*TAN (EPS) / (L* (H+S4-—-L3) ) 


TF (CHHL.EQ.2 .OR. CHSB.EQ.2 .OR. CHSB.EQ.4) THEN 


oF. 


CY=(S2*L3*TAN (EPS) ) / ((S2—-L1) * (H+S4-L3) ) 
ENDIF 


DY=S2*TAN (EPS) / (S*K) 

EY=0 

FY=0 

AZ=L5/ (S2+L1) +L3/ (S*TAN (EPS) ) 
BZ=S2*L5/ (S*K) 
CZ=(L2/L) * (L3/ (H+S4—-L3) +1) 


IF (CHHL.EQ.2 .OR. CHSB.EQ.2 .OR. CHSB.EQ.4) THEN 
CZ=(S2/ (S2-L1) ) * (L3/ (H+S4—-L3) +1) 
ENDIF 


DZ=S2/ (S*K) 
EZ=1/ (S2+L1) 
FZ=1/ (S*TAN (EPS) ) 
LXFWD=P XC 
LYFWD=BY*PYC+CY*PZC+DY *MXC 
LZFWD=AZ *PXC+BZ*PYC+CZ*PZC+DZ *MXC+EZ *MYC+FZ*MZC 
LTOTFWD= (LXFWD**2+LYFWD**2) **.5 
WRITE (NUMOUT, *) 
WRITE (NUMOUT,*)’ FORWARD HOOK/LUG LOADS LBS’ 
WRITE (NUMOUT, *) ‘ VERTICAL =! , LZFWD 
WRITE (NUMOUT, *) ’ LATERAL =’, LYFWD 
WRITE (NUMOUT, *)’ LONGITUDINAL =’, LXFWD 
WRITE (NUMOUT, *) ’ TOTAL SHEAR =', LTOTFWD 
WRITE (NUMOUT, *) 
NUMOUT=6 
WRITE (NUMOUT, *) 
WRITE (NUMOUT,*)’ FORWARD HOOK/LUG LOADS LBS’ 
WRITE (NUMOUT, *) ’ VERTICAL =’, LZFWD 


WRITE (NUMOUT, *) ’ LATERAL =! , LYFWD 
WRITE (NUMOUT, *) LONGITUDINAL =’, LXFWD 
WRITE (NUMOUT, *) TOTAL SHEAR =', LTOTFWD 

NUMOUT=11 

RETURN 

END 


C* SUBROTINE TO CALCULATE AFT HOOK LOAD 
SUBROUTINE AHOOK 
INCLUDE ’TOPJUNK.INC’ 
AY=0 
BY=(S1/S) * ((L5* (TAN(EPS) ) /K)-1) 
C* FIND SIGN OF BY 
I=S5*TAN (EPS) — (S3+L3) 
IF (I. LT.0) THEN 
BY=-BY 
ELSEIF (I.GT.0) THEN 
BY=BY 
ELSE 
BY=0 
ENDIF 
CY=L1*L3*TAN (EPS) / (L* (H+S4-L3) ) 


IF (CHHL.EQ.3 .OR. CHSB.EQ.2 .OR. CHSB.EQ.7) THEN 
CY=(S1*L3*TAN (EPS) ) / ((S1-L2) * (H+S4-L3) ) 
ENDIF 

DY=S1*TAN (EPS) / (S*K) 

EY=0 

FY=0 
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AZ=L5/ (S$1+L2)+L3/ (S*TAN (EPS) ) 
BZ=S1*LS/ (S*K) 
CZ=(L1/L) * (L3/ (H+S4-L3) +1) 


IF (CHHL.EQ.3 .OR. CHSB.EQ.2 .OR. CHSB.EQ.7) THEN 
CZ=(S1/ (S1-L2) ) * (L3/ (H+S4—-L3) +1) 

ENDIF 

DZ=S1/ (S*K) 

EZ=1/ (S1+L2) 

FZ=1/ (S*TAN (EBS) ) 

LXAFT=PXC 

LYAFT=BY*PYC+CY*PZC+DY*MXC 

LZAF T=AZ *PXC+BZ*PYC+CZ*PZC+DZ *MXC+EZ *MYC+FZ*MZC 
LTOTAFT= (LXAFT**2+LYAFT**2) **.5 

WRITE (NUMOUT, *) 

WRITE (NUMOUT,*)’ AFT HOOK/LUG LOADS LBS’ 


WRITE (NUMOUT, *) ‘ VERTICAL =", LZAFT 
WRITE (NUMOUT, *) ‘ LATERAL =! , LYAFT 
WRITE (NUMOUT, *) ‘LONGITUDINAL =’, LXAFT 
WRITE (NUMOUT, *)’ TOTAL SHEAR =", LTOTAFT 
WRITE (NUMOUT, *) 
NUMOUT=6 
WRITE (NUMOUT, *) 
WRITE (NUMOUT, *)’ AFT HOOK/LUG LOADS LBS’ 
WRITE (NUMOUT, *) ‘ VERTICAL =! , LZAFT 
WRITE (NUMOUT, *) ‘ LATERAL =', LYAFT 
WRITE (NUMOUT, *) ’LONGITUDINAL =’, LXAFT 
WRITE (NUMOUT, *)’ TOTAL SHEAR =', LTOTAFT 
NUMOUT=11 
RETURN 
END 
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